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Abstract: We discuss the possibility that flavor hierarchies arise from the electroweak 
scale in a two Higgs doublet model, in which the two Higgs doublets jointly act as the 
flavon. Quark masses and mixing angles are explained by effective Yukawa couplings, 
generated by higher dimensional operators involving quarks and Higgs doublets. Modihed 
Higgs couplings yield important effects on the production cross sections and decay rates of 
the light Standard Model like Higgs. In addition, flavor changing neutral currents arise at 
tree-level and lead to strong constraints from meson-antimeson mixing. Remarkably, flavor 
constraints turn out to prefer a region in parameter space that is in excellent agreement 
with the one preferred by recent Higgs precision measurements at the Large Hadron Collider 
(LHC). Direct searches for extra scalars at the LHC lead to further constraints. Precise 
predictions for the production and decay modes of the additional Higgs bosons are derived, 
and we present benchmark scenarios for searches at the LHC Run IL Flavor breaking at 
the electroweak scale as well as strong coupling effects demand a UV completion at the 
scale of a few TeV, possibly within the reach of the LHC. 
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1 Introduction 

The origin of the observed hierarchies in fermion masses and mixings remains one of the 
most intricate puzzles of the Standard Model (SM) of particle physics. The sizes of the 
Yukawa couplings range over at least six orders of magnitude, and the magnitude of the 
CKM matrix elements varies between 1 and 10“^. Various extensions of the SM have been 
proposed in order to explain these hierarchies. In a seminal paper, Froggatt and Nielsen 
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introduced an abelian flavor symmetry by which only the top Yukawa coupling is allowed 
as a renormalizable operator [1]. The remaining Yukawa couplings are generated as higher 
order effective operators, schematically given by 

0=y(f) QHqn, (1.1) 

where lighter fermion masses require additional insertions of the Froggatt-Nielsen scalar, or 
flavon S. At a given energy scale, the flavon acquires a vacuum expectation value (S) = / 
and breaks the flavor symmetry. The fundamental Yukawa couplings y are anarchic and 
hierarchies in the effective Yukawas are generated by the exponents re of the ratio //A < 1, 
where A is the scale at which new physics sets in. While the Froggatt-Nielsen paradigm 
does neither specify the flavor breaking scale / nor the new physics scale A, the later imple¬ 
mentation of this mechanism by Babu and Nandi [2] and Giudice and Lebedev [3] relate the 
flavor breaking scale to the electroweak scale. In particular, they propose S'/A —?> WH/h? 
in (1.1). This interesting idea however has the shortcoming that the bilinear H is a 
singlet under all symmetries, in particular it cannot carry a flavor charge. As a conse¬ 
quence, the number of flavon insertions needed in order to generate the observed fermion 
mass hierarchies is ad hoc and not related to a flavor symmetry. As briefly mentioned in 
[3], such a connection between the electroweak and the flavor breaking scale can however 
be motivated in a supersymmetric model featuring two Higgs doublets. Phenomenological 
constraints from the SM Higgs mass and signal strengths measurements exclude both the 
original Babu-Nandi-Giudice-Lebedev model as well as a possible (minimal) supersymmet¬ 
ric extension. 


A 

/ 


A 

u w / 


Figure 1: Mass scales in a generic Froggatt-Nielsen model (left) compared to the model 
proposed here (right). 

In this article, we propose a two Higgs doublet model, in which the two scalars 
and Hd act jointly as the flavon field, such that S/A —)• As a consequence, the 

flavor breaking scale is set by the electroweak scale, u ~ /, and the new physics scale is in 
the ballpark of a few TeV, as sketched in Figure 1. 

In the present study we concentrate on the quark sector and include the tau Yukawa 
couplings, reserving a full treatment of the lepton sector for future work. We discuss 
Higgs phenomenology, as well as its connection to flavor physics and show the potential for 
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distinctive discovery signals that point towards an explanation of flavor at the electroweak 
scale. 

In our model, the Higgs dependent effective Yukawa couplings induce tree-level flavor 
changing neutral currents (FCNCs) mediated by the Higgs bosons. These FCNCs, although 
naively very large, turn out to be under control for a sizable region of the parameter space. 
To this end we perform a careful study of FCNC effects in K — K, ~ g mixing 
and estimate effects in the inclusive Bg —)• Xg'y decay as well as in the flavor-violating 
top decay t —>■ he. Flavor diagonal couplings of the SM-like Higgs to quarks, as well 
as couplings between the Higgs and electroweak gauge bosons, are modified with respect 
to the SM. While the former are unique to our model, the latter are equivalent to the 
Higgs couplings to gauge bosons in generic two Higgs doublet models [4, 5]. This leads 
to deviations in both the Higgs production cross section and decay rates and we compute 
these effects for all relevant channels to compare them with current bounds from both 
the ATLAS and CMS experiments. We perform a global fit to all SM Higgs LHC data 
and we can accommodate the experimental data at a 2a level for a sizable range of model 
parameters. It is most remarkable, that the parameter space preferred by flavor observables 
has a signiheant overlap with the region preferred by the SM-like Higgs global fit. 

A characteristic feature of this two Higgs doublet flavor model is, that both the con¬ 
straints from Higgs signal strength measurements and flavor physics point to a parameter 
region far from the alignment/decoupling limit, such that the additional Higgs bosons 
cannot be arbitrarily heavy. Furthermore, electroweak precision observables favor a large 
mass splitting between charged and neutral scalars, while the neutral scalar masses are 
preferred to be almost degenerate. As a result, direct collider searches for the additional 
Higgs bosons are very powerful in probing this model. We analyze the LHC results from 
direct searches for the CP-even and CP-odd Higgs scalars as well as for the charged Higgs 
boson in various production and decay modes and identify the most promising channels for 
a discovery. Although the bosonic Higgs couplings parametrically correspond to the ones 
in a generic two Higgs doublet model, the parameter space singled out by flavor constraints 
and Higgs precision measurements leads to distinctive predictions for future searches at the 
LHC. 

Altogether, the two Higgs doublet flavor model presented in this work provides an 
explanation for quark masses and mixing angles from physics at the electroweak scale, 
while providing new opportunities for Higgs phenomenology at the LHC. The model can 
be tested by high precision measurements of meson-antimeson mixing and implies a UV 
completion at a scale that can be probed at the LHC. 

This paper is organized as follows: In Section 2 we introduce our model, discuss the 
relevant parameters in the Yukawa sector and constraints from quark masses and mixing 
angles. We subsequently compute the Higgs couplings to quarks in Section 3. In Section 4, 
5 and 6 we investigate constraints from Higgs, flavor and electroweak precision observables 
and map out the parameter space in agreement with these constraints. Section 7 contains a 
detailed analysis of present and future collider searches for the extra scalars. We comment 
on a possible UV completion Section 8. In Section 9 we present benchmarks for our model, 
before we summarize our main results in Section 10. 
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2 Flavor from the Electroweak Scale 


We consider a two Higgs doublet model in which fermion masses are generated by a 
Froggatt-Nielsen mechanism. We assume that the combination of the two scalar doublets 
HuHfi carries a non-zero flavor charge such that the flavon is replaced by 

A2 - a2 • 

We assign opposite hypercharges to the two Higgs doublets and parametrize them as 


Hu = ^ 


(vu + ReH^ + i\iaHl 

\/2 V V2H- 




Vd + Re + i Im 


( 2 . 2 ) 


In this setup the electroweak scale sets the flavor breaking scale by 

/ ^ {HuHd) _ VuVd 

A A 2 2 A 2 ’ ^ ^ ^ 

where 

= + —= tan/I, (2.4) 

Vd 

with V = 246 GeV and 0 < /3 < 7 r/ 2 , such that Vu and Vd are positive. We define the 
expansion parameter 

^VuVd ^ tan/3 

2A2 “ l+tan^2 2A2 ■ 

We choose e = rrih/mt ~ 1/60, such that the Yukawa coupling for the bottom quarks 
corresponds to an effective operator with one insertion of the Higgs doublets (n = 1 in 
terms of equation (1.1)). Therefore for tan/3 = 1, the new physics scale is approximately 
A ~ 4n ~ 1 TeV. If the fundamental Yukawa couplings in the UV completion are slightly 
larger than 1, this bound becomes weaker, and values of tan/3 > 1 are possible with a UV 
scale of the order of a TeV. Therefore, an ultraviolet completion at the TeV scale and tan /3 
of 0(1} are predictions of this model. We further discuss such a UV completion in Section 8 . 


We consider the quarks and scalars in our model to be charged under a global U{1 )f 
symmetry. Therefore in the flavor eigenbasis the Yukawa sector of the SM is replaced by 
the effective Lagrangian (to leading order in powers of e) 


-^Yuk = y, 




HuH, 

A2 


Qi^uVjdj 4 “ 


V 


HuH, 

A2 


QiHdduj + h.c ., ( 2 . 6 ) 


in which = au,ac,at, and ad^ = ad, ag, at denote the flavor charges of the three gener¬ 
ations of up- and down-type quark singlets, a* = 01 , 02,03 the flavor charges of the three 
generations of quark doublets and oh„, the flavor charges of the Higgs doublets. The 
leading order Yukawa couplings in equation (2.6) reduce to the Yukawa sector of a two 
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Higgs doublet model of type II in the limit of vanishing flavor charge Oj, Uuj , adj —)• 0. Cou¬ 
plings of Hu{Hd) to the down- (up-) type quarks are suppressed by additional powers of 
The fundamental Yukawa couplings yf^ and yf^ are considered to be anarchic and of 0(1). 
In writing equation (2.6) we normalized the sum of the Higgs charges to + clh^ = 1. 
The effective Yukawa couplings are then given by 


iyu)ij = ytj , iYdhj = yfj • (2.7) 

In (2.6) and (2.7), repeated indices between ytj and are not summed over, i.e., 

for example (I7i)i2 = ■ Thus the hierarchy of the effective Yukawa couplings 

is determined by the structure of the exponents of e. The rotation to the mass eigenbasis 
is performed via 

yu,d = Uu4 K,d ’ (^•^) 


with diagonal matrices given by 


A 


U 


V2 

Vu 


diag{mu,mc,mt) , 


Ad 


Vd 


diag{md,ms,mb) , 


(2.9) 


and unitary rotation matrices Uu^d, Wu,d- 

In the following we fix the flavor charges of the quarks and Higgs bosons by imposing 
constraints from quark masses and the CKM matrix. If the charges of the three generations 
of quark doublets and singlets are ordered such that 


ttl ^ 0 j 2 ^ U3 5 Y (Xd , 


( 2 . 10 ) 


one can derive the 0(e) dependence for the quark masses and rotation matrices [1], 


'^u a, — Ou . —a tf,, 

niuj OC —= £ ^ 5 

V 2 


'^d aj—ad—aij, 

rudj (X —^e ^ “d 

v 2 


(t/,)i,-(xel““-“^'l, (IT„)i,-(xel““*-““^'l, (ITd).^ oc 


( 2 . 11 ) 


for i,j = 1,2,3. In the numerical analysis we will use the full unitary rotation matrices 
and include a scanning of anarchic Yukawa couplings with arbitrary phases and absolute 
values \yy^\ G [0.5,1.5]. Six of the 11 flavor charges are fixed by the quark masses. We 
choose 


Vu 

mb ^ 

^ nic ^ 

1 

ms ^ 

2 

rud ^ 



mt 

mt 

- ^ , 

mt 

, 

mt 

mt 


Additional conditions follow from the CKM matrix. 


( 2 . 12 ) 


VcKM = Ul^Ud, (2.13) 

^We also explored choices of flavor charges in which both up- and down-type quarks couple to one of the 
Higgs doublets at leading order (based on a two Higgs doublet model of type I), which will be discussed in 
a separate publication [6]. 
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by fixing 


(Vckm)i2 ~ , (Vckm)i3 ~ (Vckm)23 ~ • (2-14) 

These conditions end up fixing only two parameters. Including the normalization of the 
Higgs charges + clh^ = 1 and our choice of = 1, we have 10 conditions on 11 
parameters^. The remaining choice allows for an overall shift of quark flavor charges. 
Physical quantities however only depend on invariant differences. Thus the remaining 
choice does not have any phenomenological consequences and we set 




Oi 

= 2, 

— 

-2, 

ad — 

-1 

02 

= 2, 

dc — 

0, 

as — 

0 

03 

= 1, 

at = 

0, 

ab — 

0 


If the last condition (2.14) is replaced by 

(Vckm)i2 ~ (VcKm)i3 ~ (Vckm)23 ~ , 


(2.15) 


(2.16) 


only the structure of the quark masses is explained by the flavor charges, while the hier¬ 
archical form of the CKM matrix is determined by the fundamental Yukawas yfj. In 
this case, a suitable choice of flavor charges read 


= 1 , 

aHa = 0 , 


Ol 

= 2, 

au — 

-2, 

O-d — 

-1 

02 

= 2, 

ttc ~ 

0, 

as = 

0 

03 

= 2, 

at = 

1, 

Ofo 

1 


(2.17) 


This choice of charges is motivated by considerably weaker constraints from flavor observ¬ 
ables due to the aligned charges for the left-handed quark fields. 


A detailed implementation of lepton masses and mixing angles is beyond the scope of 
this work. We will however define the couplings of the tau leptons to the scalars in our 
model, since they are important for the Higgs phenomenology. We set 

0^, = yr^^^fLHdTR, (2.18) 

such that mr/mt ~ e. 


3 Higgs Couplings 


The Yukawa interactions give rise to modifications to flavor diagonal Higgs couplings as 
well as potentially dangerous flavor changing neutral currents. In the flavor eigenbasis the 
interaction between quarks and the real neutral components of the Higgs doublet scalars 
follows from (2.6) and we obtain 


>Co = {yu)ij [(1 + Oi - Onj - ana) Rei7° -|- (a, - - ano) tan /3 Re 17^] ul,ur. (3.1) 


(Yrf)-j. (l + Ui- - a ho) Re + (a* - ad^ - Oho) cot /3 Re H) 


d 


d-LidR. +h.c.. 


^Different choices for the normalization condition or the Higgs charges, e.g. an^ = = 0, do not 

change the physics of this model but will only imply different assignments for the quark flavor charges. 
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We rotate to the quark mass eigenbasis, according to equation (2.8) and introduce the 
Higgs mass eigenstates as defined in Appendix A. The rotation of the scalars gives rise to 
the following couplings between the scalar mass eigenstates and quarks 


A = (bu) ULi h UR^ + {hd) dii h H ur. + (Bd) di, H dR^ + h.c., (3.2) 

in which 

{bu)ij = (Tu).^. [(1 + Oi - - anjcosa- (a* - - anj sinatan^] , 

{bd)ij = {yd)ij [-(1 + Oi - adj - ttHa) sin a + (oj - adj - auj) cos a cot /?] , 

{Bu)ij = [(1 + Oi - - aH„) sina + (a* - - anj cosatan^] , 

{Bd)ij = iYd)ij [(1 + Oj - ttdj - aui) cos a + (a, - adj - anj sin a cot fd] . (3.3) 

After rotating to the quark mass eigenbasis, 


dhuiUj — {u}j)ik {bu)kl {^u)lj 1 dHuiUj — {Uu)ik {Bu)kl i^u)lj ; 

9hd,dj = {Ud)ik {bd)ki {y^d)ij , QHd^dj = {ul)ik {Bd)ki {Wd)ij , (3.4) 


we find for the couplings of the light neutral scalar h, 


Qhdidj — 


/ruu^ 

\ r 

Iet 

).. "b 


ij 

frrid^ 

1 A 


1.. bij 

\ V / 

V 


— - a.H^ f{a,f3) 


+ /(«,/3) 


—- -o,Ha /(a, fd) 
C/3 


and for the heavy neutral scalar H 

s, 


dHuiUj — I 


dHdidi — 


/ruu' 

1_ 

to 

V V 


(rud'' 

1 A 

Itt 

1 Oij 

\ U / 



- F{a,ld) 


y^-aH,F{a,fd) 


+ f (a, fd) 

+ Fia,(d) 
+ F{a,(d) 


■v \ 

Q: 




m„ 


d 


V / a 

i n 


^ V / jj \ V / a 






/rrhu 

V nj 


rriy 


U, 


=?.(?),-(v) 


V / a 

i n 


V 




(3.5) 


(3.6) 


in which = diag(mtj, tUc, m^), rud = diag(m(i, ruf,) and we define Stp = simp^c^p = 
cost/? and t^p = tany?, for any angle ip. In both (3.5) and (3.6), repeated indices are not 
summed over and we suppress the chirality index of the fermions qi = qr^, qj = qR^. We 
make use of the following trigonometric functions 


f{a, Id) = — - — = Cp-a (r , 

S/3 cp \ti3 ) 

F{a, jd) = - 1 -= 2c^_„ + s^_Q, (tp — —^ , 

C/3 S/3 H M 


(3.7) 


which are universal for up- and down-type quarks. We also define the matrices 


£=l £=l 

3 3 

^b = E = E 


(3.8) 


k=l 


k=l 
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The structure of these matrices is fixed by the flavor charges, as given at the end of Section 
2. We find for the flavor charges in (2.15), 






(2 £2 A 


(-2 £2 £ 2 \ 


^—1 £ £\ 

£2 2 £ 


£2 £2 £4 


£ £2 £2 

\ £ £ ij 


1^£2 £4 e^j 


(^ £ £2 £ 2 y 


(3.9) 


For completeness, we also give the expressions for these matrices in the case of the flavor 
charges (2.17), 


to 

o 

o 


(-2 £2 £3\ 


1 £ £2\ 

0 2 0 


£2 £2 £ 


£ £2 £ 

^ 0 0 2y 


^ £3 £ l) 


1^£2 £ ij 


(3.10) 


Note that all flavor off-diagonal Higgs couplings are proportional to these matrices. In the 
limit of degenerate flavor charges a*, am or ad^, these matrices become diagonal and do not 
induce any flavor violating couplings. For the flavor charges (2.17), therefore only U and 
V generate FCNCs. 

In addition, all flavor violating couplings of the scalars in (3.5) and (3.6) are propor¬ 
tional to the trigonometric functions in (3.7). In the limit /(a, /?) = 0, all flavor off-diagonal 
couplings of the light Higgs vanish and the diagonal couplings are independent of both CjS-a 
and tjs, and approach their SM values (up to a sign). It should be noted that this sign 
difference corresponds to the wrong-sign Yukawa coupling in a generic two Higgs doublet 
model [7, 8]. We will come back to these observations when we discuss flavor observables 
in Section 5. The limit c^_o = 0, associated with decoupling [9-11] or alignment [5, 10-12] 
is not the SM, but corresponds to the model proposed by Babu, Nandi [2], and Giudice 
and Lebedev [3]. 


The pseudoscalar mass eigenstate A is obtained through the rotation (A.3) and its 
couplings to quark mass eigenstates can be derived from (3.6), by replacing 


dAqiQj 


i dHqiQj 


CoL—^Sj^^ Sa—^Cp 


(3.11) 


Finally, the charged Higgs couplings can also be obtained from (2.6) and are indepen¬ 
dent of the flavor charges. After rotation to quark and Higgs mass eigenstates, see (A.4), 
we obtain 


r V 2 1 . 

C± = —— (m^ 

V 


Ikj 


Vr 


t 


CKM 


ik 


y/2 

dii H~ UR. H - {uud)kj (^CKM)jfc ULi dR. + h.c. 


(3.12) 


The couplings of the charged Higgs to quarks are therefore equivalent to the ones in the 
two Higgs doublet model of type H, see for example [13]. 
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4 Higgs Production and Decay 


A light SM-like Higgs has been discovered at the LHC in various decay channels. While 
observations are mainly in the ballpark of SM expectations, there is still room for new 
physics. The modihed flavor diagonal fermion couplings of the light Higgs h introduced in 
the previous section as well as modihed gauge boson couplings lead to deviations in both 
production cross section and decay rates. In the following we compute these deviations and 
compare the results with the proton-proton collision data at ^/s = 7 and 8 TeV obtained 
from the ATLAS [14] and CMS [15] experiments. 


For a given Higgs boson production channel and decay rate into specihc hnal states 
X, normalized to the SM values, we dehne the signal strength parameter 




<^prod 

^SM pSM 
"prod h^X 


■pSM 
^ h, tot 


Th 


(4.1) 


New physics can enter each of these three quantities: the production cross section iTprod, 
the partial decay rate Th^x and the total width T/j^tot- We quantify the changes in flavor 
diagonal couplings of the Higgs to fermions / = t,5, r and to vector bosons V = W^,Z 
with respect to the SM by 


SM 

9hff = HfQhff = ^f — 


ghvv = i^v 9hvv = 


2my 


(4.2) 


such that Kf = Ky = 1 in the SM limit. 

It follows from equation (3.5), that the coupling of the light Higgs to the top quark is 
rescaled by 


Cq Cfj—a 

Kt = — = — -h Sjs-a . 


(4.3) 


As a result, these couplings are modihed in the same way as in two Higgs doublet models of 
type H, see for example [5, 13, 16]. However, couplings to the other flavors signihcantly dif¬ 
fer from the couplings in generic two Higgs doublet models because of the Higgs dependent 
effective Yukawas, such that the Higgs-bottom coupling is rescaled by 


Kb — —2 - 1 -— "isp-a + C/3_a f - - . (4.4) 

Note, that for f{a, /3) = 0, any dependence on cp-a and tp cancels in (4.3) and (4.4) and 
we find that Kt = 1 and Kb = —1 and therefore the light Higgs has couplings to fermions 
of SM strength. We illustrate the parameter dependence of the square of these couplings 
in Figure 2. In the right panel of Figure 2 the value of goes through zero signalizing Kb 
changes sign and becomes negative in the upper right (lower left) corner for cos(/3 — a) > 0 
(cos(/3 — a) < 0). The structure of these couplings has signihcant impact on the Higgs 
boson production cross sections and decay rates. Further, the coupling of the light Higgs 
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Cos(j6-a) Cos(j6-a) 

Figure 2: Contours of nj (left) and (right) in the cos(;0 —a) — tan/3 plane, nf = = 1 

corresponds to the SM limit, for m, up to a sign in the right upper (lower left) corner 
for cos(/3 — a) > 0 (cos(/3 — a) < 0). The decoupling/alignment limit corresponds to the 
Babu-Nandi-Giudice-Lebedev model. 


boson to charm quarks is rescaled by 

Kc = + Cp-a . (4.5) 

In general, fermion mixing effects generate corrections to the couplings, since the flavor 
charges of the quarks are not universal. These effects are encoded in the matrices lA 
and T) given in equation (3.8). For flavor-diagonal Higgs couplings to fermions we neglect 
corrections of O (e). For couplings of the light Higgs boson to tau leptons we assume that a 
mechanism similar to our findings in the quark sector is responsible for generating masses, 
such that 


Kr = Kb . (4.6) 

For the couplings of the light Higgs to vector bosons we obtain 

Kv = , (4.7) 

which is the same as in generic two Higgs doublet models. 


The gluon fusion initiated Higgs production, neglecting light quark contributions in 
the fermion loops, is defined normalized to the SM value as 


^gg^h 
^gg^h 


l + Cfe 


Kb 

Kt 


(4.8) 
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Decay Mode 

Production Channels 

Production Channels 

C^VBF, crvH 

Experiment 

h^WW* 

fiw = [17] 

HW ^ 0.75 + 0.35 [18] 

Hw = 1-2710;®3 [17] 

Hw - 0.7 + 0.85 [18] 

ATLAS 

CMS 

h ZZ* 

^^z = 1.7l°i [19] 

= 0.81°:^® [20] 

= O.Stli [19] 

HZ = 1.7+1:! [20] 

ATLAS 

CMS 

h ^ 'j'j 

= 1.32 + 0.38 [21] 

= 1.13l°j[ [22] 

Hj = 0.8+ 0.7 [21] 

[22] 

ATLAS 

CMS 

h^bb 

^lb = 1.5 + 1.1 [23] 

+6 = 0.67l};i [25] 

/r6 = 0.52 + 0.32+ 0.24 [24] 

Hb = 1.0+ 0.5 [26] 

ATLAS 

CMS 

h ^ TT 

Hr = 2.0 + 0.8+o;8 + 0.3 [27] 
Hr ^ 0.5t[];8 [28] 

Hr = l.2At^oil -o:29±0.08 [27] 
[28] 

ATLAS 

CMS 


Table 1: Input data for the global of Higgs production and decay with references. 

The data includes all updated results of the pp collision data at y/s = 7 and 8 TeV obtained 
from the ATLAS [14] and CMS [15] experiments. 


where = —0.032+ 0.035 i depends on the loop functions given in [4]. Therefore for values 
of Kfc of 0(1), the main Higgs production channel is to leading order indistinguishable from 
a type H two Higgs doublet model. Vector Boson Fusion (VBF) and Higgsstrahlung (VH) 
are both rescaled by Ky, while associated Higgs boson production with a top pair is modified 
by Kt, 


^tt^h _ ^2 

_S_M “ 


and 


O'VBF 


fjSM (T 
OvBF o- 


O'VH _ 2 

SM — '^V ■ 


VH 


(4.9) 


Therefore the three production processes rescale with the same factors as in generic two 
Higgs doublet models, as given e.g. in [5, 13, 16]. 

The partial decay widths of the light Higgs into SM fermions / and gauge bosons 
V = , Z can similarly be written as 


j 

t^SM 


'^h^VV _ ,J2. 
pSM “ • 

^ h^VV 


(4.10) 


Both top quark and boson loops enter the diphoton decay width [29], 


pSM 

h^'y'y 


0 . 28 Kt — 1 . 28 ku / + S 


(4.11) 


in which contributions from light fermions are neglected and contributions from charged 
scalar loops are encoded in 5. We find for > 300 GeV a contribution of less than 

6 < 0.04 and set it to zero in the following [9, 29]. 
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Figure 3: Allowed Icr (dark red) and 2a (light red) regions, for a global fit to ATLAS 
and CMS data from measurements of Higgs boson decays in the left and right panel, 
respectively. The channels entering the fit are summarized in Table 1 and errors are 
symmetrized. 


Expressed in terms of the rescaling factors Kt, Kb, Kc, Kr and Ky, the total Higgs boson 
width is given by [30, 31] 

p 

« 0.57 kI + 0.25 + 0.09 k^ + 0.06 + 0.026 k^ + 0.004 , (4.12) 

where T®^ = 4.07 MeV [32] and we assume h —)• Zj and even rarer modes to be SM-like. 
These contributions are collected in the constant term 0.004. 

The partial decay width into bottom quarks has a very different dependence on tan /3 
and cos(/3 — a) than in the generic type H two Higgs doublet model. This plays a relevant 
role in dehning the allowed region in parameter space, since the bottom quark partial decay 
width dominates the total decay width, that in turn importantly affects the signal strength 
for all channels. 

In Figure 3 we show the result of a global fit based on the data collected in Table 
1. Symmetrized errors are used for the ht. The left panel shows the plot for ATLAS and 
the right panel the plot for CMS. The two fit parameters are cp-a and tp. The Icr and 
2a regions consistent with the LHC data are shaded in dark and light red, respectively. 
It is clear, that the preferred parameter space is different from generic two Higgs doublet 
models, for which regions close to the alignment or decoupling limit c^_q, = 0 are favor¬ 
able. [5, 11, 33]. In our case, c^_a = 0 corresponds to the Babu-Nandi-Giudice-Lebedev 
model [2, 3], which is clearly disfavored by the data. We observe, that while the allowed 
region for ATLAS is slightly smaller than in the case of CMS, both hts show a preference 
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Production: VBF, VH 



Cos(P-a) 


Cos(j6-a) 


Figure 4: The upper panels show the signal strengths fj,]y (left) and fiz (right) and the 
lower panel the signal strengths (left) and fXb (right) plotted against c^-a- The red (blue) 
band is the symmetrized Icr region of the corresponding ATLAS (CMS) measurement. Each 
plot shows curves for = 3 (solid red), t/j = 2 (dashed orange), tjs = 1 (dot-dashed green) 
and tp = 0.5 (dotted blue). 


for values of c^_a > 0 and > 1. The more constrained region of parameter space for 
ATLAS can be understood by the larger central values of /r^, and /r.y in the dominant 
gluon fusion channel, that are less compatible with larger values of Kb, see Figure 2. The 
white area between the two branches in both fits can be explained by very small values 
of Kb for which all other branching fractions grow. Overall, the fermion couplings prefer a 
region in parameter space, where they approach their SM values, with the caveat that the 
value of the bottom Higgs coupling Kb has a negative sign with respect to the SM value 
in the upper right branch of the allowed red region. Note also that small values of c^_o 
correspond to larger tp in the region preferred by the global ht as follows from equation 
(4.4). 

In order to understand the features of the global fit, we present the signal strengths of 
the relevant decay channels in Figure 4. In these plots, the red (blue) band is the la region 
of the corresponding ATLAS (CMS) measurement. Each plot shows the prediction of a 
particular signal strength for /rvV) Mz, fJ^-y and fib, depending on CjS-a for tjs = 3 (solid red), 
t /3 = 2 (dashed orange), = 1 (dot-dashed green) and tjs = 0.5 (dotted blue). Excluding 
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Figure 5: Tree-level contributions to AS* = 2 processes. 


all but these four observables only marginally changes the global fits. For > 1 all four 
measurements prefer values of c^-a > 0. There is also an allowed region for c^_a < 0 
for values of < 1, however as will be shown later this region is phenomenologically less 
interesting. 

We conclude, that the global fit to LHC Higgs measurements accommodates tan/3 of 
0(1) for sizable values of cos(/3 — a) away from the decoupling/alignment limit. This is a 
nontrivial result, given that tan /3 is already constrained to be of order one from the bound 
on the new physics scale. As we discuss below, values of tan /3 < 5 are in agreement with 
flavor constraints as well as a possible UV completion scale in the TeV to a few TeV range. 


5 Constraints from Flavor Observables 


In addition to modifications of flavor-diagonal couplings, the misalignment of the mass 
and coupling matrices induces flavor changing couplings of the light Higgs h, the heavy 
neutral scalar H and the pseudoscalar A. These couplings generate FCNCs at tree-level, 
which are subject to strong constraints from neutral meson oscillations. In the following 
we calculate and analyze contributions to the relevant observables. We further estimate 
effects in 6 —)■ S 7 and give the prediction for the flavor-violating top decay t —)■ he. 


5.1 Meson-Antimeson Mixing 

In the K — K system, contributions from Higgs mediated FCNCs are captured by the 
effective Hamiltonian 


^NP ^ ~ (sl 7/1 dL)^ + (sr 7 ^ (sr dL)‘^ + (sl dR)"^ 

+ Cf (sr dL) (sl dR) -h (sl 7/i dL) (sr 7^dR) + h.c. . (5.1) 


At tree-level, the corresponding Wilson coefficients can be read off from the diagrams in 
Figure 5 [34], 
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(5.2) 
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Table 2: Flavor specific part of the Wilson coefficients for meson-antimeson mixing in 
the case of the flavor charge assignments (2.15) with flavor structure (3.9) (left) and flavor 
charge assignments (2.17) with flavor structure (3.10) (right) . 


Similar expressions hold for Bg — Bg mixing, with sd —)• 6s, B^ — B^ mixing, with 
sd —)• bd and D — D mixing with sd uc. Contributions from Higgs boson exchange are 
only suppressed by the weak scale, but the Froggatt-Nielsen mechanism induces an addi¬ 
tional suppression of flavor off-diagonal couplings by the masses of the involved quarks as 
well as the expansion parameter e. The relative size of the Wilson coefficients (5.2) depends 
therefore strongly on the explicit flavor structure. For the flavor charge assignment (2.15), 
which is tailored to explain quark masses as well as CKM mixing angles, we collect the 
results in the left hand side of Table 2. 


In the case of K — K mixing, we find that the largest coefficient is C 2 with 
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ml 


Ml 


t0j Mlj 
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(5.3) 


where we factored out the light Higgs mass in the second line, the trigonometric functions 
/(a, (5) and F{a, j5) are defined in (3.7), and is the flavor-dependent part of the Wilson 
coefficient given in Table 2. The same expression holds for the Wilson coefficient with 
the additional suppression due to the replacement of c?^ —)• c^. The flavor-dependent 
Wilson coefficient cf' is also suppressed by e with respect to but the minus sign in the 
last line of (5.3) is replaced by a plus, which corresponds to a constructive interference of 
the different contributions. 
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(5.4) 


The limit of exact cancellation in and and maximal interference in corre¬ 
sponds to the SU (2 )L symmetric limit, in which operators of the type {siduY forbidden 
[35]. In Table 3, we present the current bounds on the Wilson coefficients at the electroweak 
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Table 3: Model-independent bounds on Wilson coefficients for meson-antimeson mixing 
evaluated at the electroweak scale in units of GeV“^ [36], taking into account the running 
described in Appendix D. The same bounds hold for the Wilson coefficients with flipped 
chirality C* —)> Q. 


scale for the different meson systems, based on [36]. These bounds have been derived by 
assuming that new physics only contributes to a single Wilson coefficient and can therefore 
only be taken as a rough upper limit. For K—K mixing, the strongest constraint comes from 
the CP violating observable ex, such that the bounds on the imaginary part of the Wilson 
coefficient is cited. Since we assume arbitrary phases, the estimate (5.3) holds for both real 
and imaginary parts of the Wilson coefficients. Comparing (5.3) with the bound in Table 
3 shows that a partial cancellation in C'|'^ is necessary in order to comply with the limit. 
For Ma-, Mh > ruh, this corresponds to a preferred region in the cos(/3 — a) — tan^d plane. 
In the left panel of Figure 7 we show the preferred region, for which lei'll < 10“^®/CeV^ 
(shaded orange), assuming Ma = Mh = 500 CeV. Contributions to C'|'^ can be enhanced 
by the constructive interference between the scalar contributions. Also, the bound on C'|'^ 
is particularly strong, because it is enhanced from Renormalization Croup (RC) running 
as well as from the matrix element, that scales like M^/(ms -|- m^)^ ~ 14, see Appendix C 
for details. However, the additional suppression shown in Table 2 gives C'|'^ = such 

that a slight enhancement from interference effects is allowed. In the left panel of Figure 
7 we show the region in the cos{/3 — a) — tan/3 plane for which < 7 x 10“^^/ CeV^ 
(shaded blue). 

In addition to tree-level exchanges, various one-loop contributions can potentially be¬ 
come large. The relevant diagrams are shown in Figure 6. The contributions from the box 
diagrams of type (a) are completely analogous to the ones in a type II two Higgs doublet 
model, because the couplings of the charged Higgs (3.12) are indistinguishable between the 
two models. The leading contribution enters and comes from the box with one charged 
Higgs [37], a boson and top quarks running in the loop and one finds 



1 1 

Ibvr^ 



1 


9 X 10”^® 


/500 CeV 

V Mh+ 


'ceV-^, 


(5.5) 


where in the last equality we set = 1. For t /3 < 1, this contribution is of the order 
of the largest tree-level contribution. We therefore require > 1 in order to be in com¬ 
pliance with experimental bounds in Table 3. In principle, there are also contributions 
from box diagrams to the other operators in (5.1) as well as box diagrams with neutral 
scalar exchange, but both are chirally suppressed by powers of light quark masses over the 
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Figure 6: Three types of one loop contributions to AS* = 2 processes. 


electroweak scale and turn out to be negligible. The loop diagrams labeled (6) and (c) in 
Figure 6 are also suppressed. Diagrams of type (6) have the same coupling structure as the 
tree-level diagrams, but are additionally suppressed by a loop factor. Diagrams of type (c) 
are enhanced with respect to (5.4) by the light Higgs couplings to the top quark or charged 
scalars, but suppressed by CKM elements and a loop factor, such that we find for C'|'^ [38] 

^4,penguin ^ 1 rnf VtsVtd ^ . ^-3 

IGtt^ e'^ f {a, P) 

The equivalent diagram with a charm quark in the loop is of the same order. 


Having considered all different contributions we will map out the parameter space in 
the cos(/? — a) — tan (3 plane in which the prediction for eK in our model agrees with the 
experimental bound within 2a in a numerical analysis. For this purpose we define 

Im(X0|^A5=2|^0)’ 

where includes the Standard Model contribution. We compute 

the Wilson coefficients at the scale of the light Higgs and for Mh = Ma = M^± = 500 
GeV respectively, using the full expressions for the Wilson coefficients including tree-level 
and leading box diagrams. We collect the full analytic expressions of the latter in Appendix 
B. In the next step, the Wilson coefficients in (5.1) are evolved down from the mass scale 
of the scalars to the scale ^ = 2 GeV at which the hadronic matrix elements are evaluated 
using the RG equations in [39]. The hadronic matrix elements are taken from [40] and 
collected with the other numerical input in Appendix D. We randomly generate a sample 
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Figure 7: The left panel shows the region in the c^_q — tg plane for which the tree- 
level contributions to lei'll < 10“^®/GeV^ (orange) and the tree-level contributions to 
< 7 X 10“^^/GeV^ (blue). In the right panel we show regions of parameter space in 
which our sample points reproduce within two standard deviations. The color coding 
indicates the percentage of points in agreement with the experimental constraint. In both 
plots, the scalar masses are Ma = Mh = M^+ = 500 GeV. 


set of points of fundamental Yukawa couplings, defined in (2.6), with G [0.5,1.5] and 
with arbitrary phases. We require the SM quark masses and Wolfenstein parameters to 
be reproduced within two standard deviations. More details to the procedure and input 
parameters can be found in Appendices C and D. At this stage, the mixing angles a and 
P from the Higgs sector still remain free parameters and our sample set only fixes the 
fundamental Yukawas. 

In the right panel of Figure 7 we show the percentage of sample points which reproduce 
within 2a in the cos(/I — a) — tan/3 plane. We employ the value extracted from a fit 
to the GKM triangle by the UTfit group [41], 

= 1.05l[];^^ @95% CL. (5.8) 

The result shows good agreement with the estimate of the separate contributions shown in 
the left panel of Figure 7. The area for which tjs < 0.5 is cut off, because of the one-loop 
contributions from charged Higgs exchange [42] . We find a large region of parameter space 
for which our model prediction is in agreement with the experimental bound without any 
tuning of parameters. 

In the case of and Bg — Bg mixing, the effective Lagrangian, as well as the 

tree-level contributions to the Wilson coefficients from scalar and pseudoscalar exchange 
can be read off from (5.1) and (5.2) with the replacements s -H- 6 and d 0 d, s, respectively. 
The angle dependence of the Wilson coefficients is universal and therefore only the flavor 
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Figure 8: In the left (right) panel we show regions of parameter space in which our sample 
points reproduce C'b^(C's^) within two standard deviations. The color coding indicates the 
percentage of points in agreement with the experimental constraint. In both plots, the 
scalar masses are Ma = Mh = M^+ = 500 GeV. 


dependent part changes from (5.3) and (5.4), such that the parametric dependence pre¬ 
sented in the left panel of Figure 7 also holds in the B sector. For the Wilson coefficients 
it follows from Table 2, 
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(5.9) 

(5.10) 


The corresponding bounds in Table 3 imply, that is at the border of the naive bound, 
while a much larger contribution to is allowed. The contributions to C^'^, and 
d'^ ciTe too large almost in the entire cos(/3 — a) — tan/3 plane, and therefore demand 
cancellations implying important restrictions for the permitted region of our parameter 
space. 

At the one-loop level, box diagrams generate the contributions 


l^bq 
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V Mh+ J 11 GeV-2 , q = s, 

(5.11) 


for tan/3 = 1. In the Bg — Bs system for low tan/3, this contribution becomes larger than 
all tree-level contributions. Since the box is only sensitive to charged Higgs couplings, we 
expect comparable constraints as in a two Higgs doublet model of type H. In addition, since 
the contribution is independent of cos(/3 — a), we expect a universal lower bound on tan/3, 
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Figure 9: Boundaries of the regions in which 10% of our parameter points agree with 
the Cb^ at the 95% CL in the positive cp-a plane. The different colors correspond to 
M = Ma = Mh = Mh+ = 400 GeV (blue) M = 500 GeV (purple), M = 600 GeV 
(green), and M = 700 GeV (light brown). 


as observed in the left panel of Figure 8. For both the — Bd^s system we also include 
the box diagram contributions to the other Wilson coefficients, which are suppressed by 
mfy/mw. The corresponding expressions are collected in Appendix B. 


Analogous to (5.7), we define 






(5.12) 


such that Csq = Amg/Am®'^ measures new physics effects in the mass difference and new 
phases enter (t>Bq- In the left (right) panel of Figure 8, we present the percentage of sample 
points in agreement with the experimental constraints at 95% CL for (C'^*^), based 
on the results obtained from the UTfit group [41], 


= 1.052l°;J^® @95% CL, @95% CL. (5.13) 


In both plots we choose Mh = Ma = Mh+ = 500 GeV. As expected from our estimate 
above, in the Bg — Bg system, we find good agreement with the experimental bounds for a 
large region of parameter space. For the Bd — Bd system, we find only a small fraction of 
the parameter space in agreement with the experimental constraints. Since the new physics 
effects in all Wilson coefficients are too large, accidental cancellations in the fundamental 
Yukawa couplings are in effect in order to achieve agreement with data. As a consequence, 
slightly tuned Yukawa couplings as well as rather heavy extra scalars Ma ~ Mh ~ 500 
GeV are necessary in order to agree with the bounds from Bd — Bd mixing. In the follow¬ 
ing, we will adopt the 10% contour as the fine-tuning bound from flavor observables on the 


- 20 







parameter space. Figure 9 shows the corresponding contours in the positive cp-a plane for 
M = Ma = Mh = Mjj+ = 400 GeV (blue) M = 500 GeV (purple), M = 600 GeV (green), 
and M = 700 GeV (light brown). The bound for low tan/3 comes from the charged Higgs 
loops in Bg — Bg mixing. A future, more precise measurement of meson-antimeson mixing 
can reveal deviations from the SM prediction or further constrain the allowed parameter 
space, if no new physics effect is found. 


hi D — D mixing, all tree-level contributions to the Wilson coefficients are strongly 
suppressed, 

TTi^ 4 x 1 n 

/^UC ^ /^UC ^ riuc '"'C ^ ^ ^ 
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(5.14) 


In contrast to the down-sector however, the box diagram with neutral Higgs exchange is 
not suppressed by light quark masses, because the dominant contribution comes from the 
top in the loop [43] . The leading box contributions of the light Higgs to the coefficient 
can therefore be larger than all tree-level effects 


~ D2{mt,mh) 

1287r2 V u (m2-m2)3 

2 X 10-16 
"" G^V^’ 


(5.15) 


for f{a,(5) = 1, and the loop function dehned in Appendix B. Boxes with heavy Higgs in¬ 
sertions are further suppressed. However, the corresponding bound in Table 3 is orders of 
magnitude weaker than our estimate. The D — D system will therefore not induce further 
constraints. 


In all the above analyses, we have concentrated on the solution for the flavor charges 
(2.15), but the situation is quite different for the flavor charges given in (2.17). From (3.10) 
it follows, that the contributions to the Wilson coefficients are highly suppressed, as is ex¬ 
plicit in the flavor-dependent parts of the Wilson coefficients given on the right hand side 
of Table 2. This shows, that although constraints from the Bg — Bg and K — K systems 
remain the same, the constraints from the B^ — Bd system can be very much relaxed due 
to the different charge assignment. Therefore, if only the hierarchies in the quark masses 
are explained by a Froggatt-Nielsen mechanism at the weak scale, but the CKM mixing 
angles have a different origin, bounds from meson-antimeson mixing are very mild and do 
not lead to any severe restrictions on the parameter space. 

Rare Kaon and Bd^g decays can in principle be subject to large corrections, but depend 
crucially on the implementation of the lepton sector, which will be discussed elsewhere. 
Processes in which the neutral scalars only enter at loop-level, such as Br{Bs —?■ Xg'y) are 
generically dominated by charged Higgs contributions, which are larger than the contribu- 
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tions from the neutral Higgs by a factor of 
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0(10^ - 10^), 


(5.16) 


for f(a,P) = 0.1 — 1. We will therefore adopt the bounds from Br(i?s —?■ Xg'j) on the 
charged scalar mass in two Higgs doublet models for tan/3 > 2, considering values within 
a 3a band in order to account for uncertainties of higher order corrections not included in 
the theoretical computation [44, 45], 


Mh± > 358 (480) GeV @ 99%(95%) CL . (5.17) 


5.2 Flavor Violating Top Decays 

We consider the flavor violating decays of the top quark t —>■ he and t —)• hu. In contrast 
to the SM, in which flavor violating top quark decays are loop suppressed, in our model 
the top quark has tree-level couplings to the light Higgs and other up-type flavors. The 
corresponding branching ratios Br(t —)• /ic) ~ 3 x 10“^^ and Br(t —>■ Lu) ~ 2 x 10“^^ are 
tiny in the SM [46]. In our model the branching fraction of the top decaying to higgs and 
charm is given by [47] 
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(5.18) 


and similarly for Br(t hu) by replacing the appropriate flavor indices. Both branching 
ratios are parametrically of the same order, because the flavor off-diagonal couplings in 
equation (3.5) yield Qhct ~ Shut oc mt£. In Figure 10 we show Br(t he) plotted against 



Cos(y6-a) 

Figure 10: The plot shows Br(t —)• he) vs. cos(/3 — a) for tan/3 = 3(4) in blue (green) as 
well as the current exclusion limits for the 8 TeV LHC (solid red) and projected limits at 
the high luminosity LHC (dashed red), respectively. 
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cos(/3 — a) for a range of parameter points and indicate the different predictions for tan /3 = 
3(4) by a blue (green) band. The widths of these bands correspond to the range of values 
obtained by scanning over our sample set of random fundamental Yukawas. The most recent 
limits are Br(t —)• he) < 0.56% from CMS [48] and Br(t —)■ he) < 0.79% from ATLAS [49] 
and are shown in the plot as a red band. The projected exclusion limit for 3000 fb~^ at the 
high luminosity LHC Br(t —)• he) < 2 x 10“^ [50] is indicated by a dashed red line. The 
plot shows that this cross section can be even above 10“^ for negative values of cos(/3 — a). 
However, the cross section drops for the same angles for which FCNCs become small, 
because the same trigonometric function governs flavor off-diagonal couplings between the 
light Higgs to up- and down-type quarks in equation (3.5). 

6 Perturbativity, Unitarity, and Electroweak Precision Measurements 

In this section we consider perturbativity bounds, as well as constraints from the unitarity 
of the S matrix and electroweak precision measurements on our model. The large scalar 
masses implied by flavor observables and the constrained scalar potential (A.l) result in 
potentially large quartic couplings. Mass splittings between the different scalar mass eigen¬ 
states can in addition generate sizable contributions to the oblique parameters 5, T and U. 
We therefore scan over the allowed parameters, considering the various bounds described 
in [51]. This includes stability constraints on the Higgs potential, perturbativity bounds on 
the quartic scalar couplings, unitarity of the various scattering amplitudes involving scalars 
and the constraints from the oblique parameters. This calculation is not different from a 
generic two Higgs doublet model, since the oblique parameters only measure corrections to 
the gauge boson self-energies from loops of the new scalars, whose couplings are fixed by 
the kinetic terms [52, 53]. 

The two plots in the upper panels of Figure 11 show the region in the positive cos(/3 — 
a) — tan /3 plane in which stability and perturbativity bounds are fulfilled, and the S and 
T parameters are at most 2cj from the best fit point, corresponding to a global fit 
obtained by the Gfitter gronp [54] . The upper left panel illustrates the allowed regions for 
degenerate scalar masses of M = Ma = Mh = Mu± = 500 GeV in light green, M = 600 
GeV in green and M = 700 GeV in dark green. For masses M = 700 GeV only valnes of 
cos(/3 — a) < 0.2 are allowed, approaching the decoupling limit. For masses M = 500 GeV 
and M = 600 GeV there is a region of parameter space in agreement with all constraints 
for values of cos(/3 — a) > 0.2, that partly overlaps the region preferred by the global fit to 
the SM Higgs signal strengths. In the upper right panel, we show the same plot for masses 
Mh+ = 360 — 700 GeV and Ma = Mh = 600 GeV {Ma = 600 GeV, Mh = 550 — 650 GeV) 
in purple (dark blue). In both upper panels, we also superimpose the 2a contours (dashed 
lines) of the global Higgs fit using the ATLAS measurements of the signal strengths, that 
are the most stringent at present. Almost all of the right branch of the global Higgs fit can 
be populated for large scalar masses, while low values of cos(/3 — a) < 0.3 are only allowed 
for tan/3 > 4.5. 

The lower left panel shows the region allowed by all constraints discussed above for 
which we further demand, that the ATLAS SM Higgs signal strengths measurements are 
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Figure 11: The upper left panel shows regions of parameter space in which the various 
constraints described in the text are fulfilled for scalar masses M = Mh = Ma = = 

500 GeV (light green), M = 600 GeV (green) and M = 700 GeV (dark green). The upper 
right panel shows the same plot for Mfj+ = 360 — 700 GeV and Ma = Mh = 600 GeV 
{Ma = 600 GeV, Mh = 550 — 650 GeV) in purple (blue). The 2a contours of the ATLAS 
fit to Higgs measurements is shown in dashed black. The lower panels show the parameter 
space in the cos(/3 — a) — Mh+ plane in agreement with all bounds discussed in the text, 
including the 2a global fit to ATLAS data. In the lower left (right) panel we assume 
Ma = Mh = 600 GeV {Mh = Ma ± (10 — 20) GeV), with values of tan j3 indicated by the 
color coding bar on the right. 


reproduced within 2a in the cos(/3 — a) — Mh+ plane for Ma = Mh = 600 GeV. The value 
of tan /3 is indicated by the color coding. The tiny gap at cos(/3 — a) Ri 0.3 is also visible in 
the upper left plot. For tan/3 < 4 only degenerate masses Ma = Mh = Mh+ or a sizable 
mass splitting of Ma — Mh+ ^ 100 GeV are allowed. We show the same plot in the lower 
right panel, but with a moderate mass splitting between the neutral Higgs boson masses. 
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Mh = Ma ± (10 — 20) GeV, while keeping Ma = 600 GeV fixed. ^ In that case the gap 
around cos(/3 — a) ~ 0.3 becomes much more prominent. 

Further, for some regions of the parameter space, one or more of the quartic couplings in 
the Higgs potential can become non-perturbative already at the TeV scale Ai(/r = ITeV) > 
Att. We implement the one-loop beta functions for our model and match to the SM at 
an approximate average scale of the Higgs boson masses in order to estimate the scale 
of strong coupling. In particular for larger values cos(/? — a) and larger and degenerate 
masses Ma = Mh, the cutoff scale becomes lower. Moreover, we hnd that for sizable mass 
splittings between the charged and neutral scalars, the scale of strong coupling is in the 
range of 2 — 5 TeV. However, as mentioned in Section 2 and in more detail in Section 8 
below, we expect the UV completion of our model to set in close to the TeV scale. 

We conclude, that for fixed Ma = 600 GeV, two qualitatively different choices of 
scalar masses are compatible with electroweak precision bounds, Higgs constraints and a 
low tan f3 as preferred by flavor constraints. Either the scalar masses are approximately 
degenerate Ma ~ Mh ^ Mh+ or the charged scalar is considerably lighter than the neu¬ 
tral scalars Ma,h — Mh+ ^ 100 GeV. Of these possibilities, only for large mass splittings 
can the theory be valid up to several TeV and in the following we will concentrate on this 
setup. Note, that these restrictions would be slightly relaxed if we take the ht to the CMS 
measurements of the Higgs signal strengths as a constraint. 


Another important electroweak precision observable is the Zbb coupling. While the 
experimental value of the left-handed Zbi^bi coupling is in good agreement with the SM 
prediction, there is a discrepancy between the measured right-handed Zbubn coupling and 
the SM prediction, see e.g. [54, 55]. Higher order corrections with the neutral or charged 
scalars in the loop can in principle affect these couplings. 

The charged scalar contributions to the left-handed Zb^bi couplings in a two Higgs 
doublet model of type H can become sizable for low tan/3, inducing a bound of tg > 0.5 
for masses of Mh± ~ 500 GeV [42] , while corrections to the ZbnbR vertex are suppressed 
by mi,/nit. In addition, the neutral scalar couplings to bottom quarks are very different 
from a generic two Higgs doublet model in a large range of parameter space. We dehne 
the couplings of the Z boson to left-handed and right-handed bottom quarks by 

J^Zbb = -75) +g^(l -75)) b, ( 6 . 1 ) 

with 

Li.R L.R I r L.R . r L.R . c L.R //^ r>\ 

9 =9sM+^9h +^9 a^h + ^9h±- ( 6 . 2 ) 

L R 

Here, are the SM couplings and we denote the corrections from neutral and charged 
Higgs exchange by 6gh, 5gA,H and 6gH±, respectively. We estimate 


^ Ml 


^9h 

-OC 

bgH± m 


TT +2 2 2 
2 ^13 ^ 

h 


^ 9 , 


R Ml -2 


5g^. 


OC 




mt 


t2 ’ 

^/3 


(6.3) 


the mass splittings become larger than \Mh — Ma\ ^ 30 GeV, the full parameter space is excluded. 
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while contributions from the heavy neutral scalars are further suppressed by 6gh/SgA,H ~ 
m\/M\ ^ and couple with and as defined in the following section in equation (7.1). 
Neutral Higgs contributions to g^ are therefore at least an order of magnitude smaller 
than the charged Higgs contributions for the region preferred by the global Higgs fit, while 
corrections to the right-handed coupling g^ are at most of a similar size. We numerically 
estimate the light neutral Higgs contributions following [56, 57]. For = 1, we find for the 
right-handed coupling Sgf^ < 10“® and for the left-handed coupling 6gf^ < 10’ 

which is many orders of magnitude too small in order to explain the anomalous Zbubji 
coupling. In order to improve the fit with respect to the SM, contributions of the order 
of 0.2% to g^yi and 2% — 20% to g^y^ (depending on the sign) are necessary [58]. The 
neutral Higgs contributions to the Zbb vertex can therefore be safely neglected. It should 
be noted, that fermionic mixing effects in the UV completion of this model can affect both 
the oblique parameters and the Zbb vertex. These however depend sensitively on the exact 
realization of the UV completion, which is beyond the scope of this paper. 

7 Collider Searches for Heavy Extra Scalars 

Our model features heavy new scalars beyond the SM, namely the neutral scalar Higgs H, 
the pseudo-scalar A and the charged Higgs H^. Their masses are bound to be less than 
700 GeV by perturbativity, and various flavor constraints set lower bounds on their masses 
as discussed in Section 5. In this section we consider the latest ATLAS and CMS bounds 
on new neutral and charged Higgs bosons. 



Cos(j6-q') 


Figure 12: Contours of in the cos(/3 — a) — tan /3 plane. A suppression of the 

coupling with respect to the SM is achieved in the darker shaded area. 
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7.1 Couplings and Total Width of Heavy Scalars 

Similar to the case of the light scalar, the couplings of the heavy scalar H and pseudoscalar 
A to quarks - with the exception of the top quark - differ from the couplings in a two 
Higgs doublet model. The couplings to gauge bosons are instead the same as in a two 
Higgs doublet model. Specifically, the couplings of H and A to gauge bosons and third 
generation quarks normalized to the SM as in (4.2), read 


H ( ^\ H 

Kfy — 3Cfj—a T Sfj—a ^ J ’ — ^13—a j 

4 = 2t^ + ^, (7.1) 

tf3 

where t, b and V denote the rescaling factor for top, bottom and vector boson couplings, 
respectively. Since is relevant for the gluon fusion production of the heavy Higgs 

boson i7, its parametric dependence is essential and we illustrate it in Figure 12. Both 
flavor diagonal and flavor changing couplings of H and A involving the charm quark, are 
given by 
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• e, 



(7.2) 


where and are defined according to equation (7.7) below. As discussed at the end 
of Section 2, we define the couplings to taus as 




= Kb ■ 


(7.3) 


The couplings of the charged Higgs 77+ to fermions are the same as in a two Higgs doublet 
model of type H. Similarly, all self-couplings between the scalars are the same as in a 
generic two Higgs doublet model. The coupling between the heavy scalar H and the light 
Higgs h is of particular interest for the following analysis and reads [5, 9, 12] 


QHhh = 


^j3—a 


(3Myi - 2mh - Mjj) I C2(/3_q,) - 


S2(J3-a) 

^2/3 


-Ml 


Finally, the couplings between two Higgs bosons and one gauge boson read [4] 
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(7.5) 


Further, we define the total widths for H, A, and 77+, including all relevant and kine¬ 
matically accessible decay channels (no off-shell decays are relevant in the regions we will 
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M = 600 GeV M = 600 GeV 



Cos(j6-a) Cos(j6-a) 


Figure 13: The plot shows the parametric dependence of the total width for the heavy 
Higgs H (left panel) and total width for the pseudoscalar A (right panel) for M = 600 
GeV. The contours, labeled in GeV, show lines of constant width. 


consider) 

rH = r{H ww) + r{H zz) + r{H hh) + t{h az) + t{h h+w-) 

+ T{H ^tt)+ T{H hh) + T{H cc) + T{H tc) + gg) 

+ T{H ^T+T-), 

Fa = F{A hZ) + r{A HZ) + r{A H+W-) + r{A ^ tt) + r{A hh) 

+ r(H —)■ cc) + r(H —)• tc) + r(H —>• gg) + r(H —)■ T~^T ~), 

Fh+= F{H+ hW+) + r(F+ ^ HW^) + r(H'+ ^ AW+) + F{H+ th) 

+ F{H+^tu). (7.6) 


Note that, besides the usual decay channels the flavor violating channel r(<l> —)• ct) with 
^ = H, A appears in 7.6. This channel is characteristic for our model and we therefore 
give the partial width explicitely 
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for M$ = Ma , 

for M$ = Mh , 
(7.7) 


with 


X{x, y, z) = x‘^ + y'^ + z'^ — 2xy — 2xz — 2yz . 


(7.8) 


The parametric dependence of the total width for the scalar (pseudoscalar) Higgs boson 
is illustrated in the left (right) panel of Figure 13 for M = Ma = Mh = Mh+ = 600 GeV. 
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M = 600 GeV 


M = 600 GeV 



Cos(j6-a) 


Cos(j6-a) 


Figure 14: We show contours of constant (j{gg —x Br(A —)• hZ) in picobarn (left 
panel) and cr(gg —)• A) x Br(74 —)■ hZ —)• iibb) in femtobarn (right panel) for 8 TeV pp 
collisions and M = 600 GeV. 


For large regions of parameter space the total width becomes large. In particular, for 
tan/3 > 1 and | cos(/3 — a)| > 0(0.5) values of 0(100) GeV can be obtained, such that 
finite width effects need to be taken into account. The charged Higgs can also have a 
sizable branching ratio Br(H'*’ —)• hW~^), which can become the dominant decay channel 
for sufficiently large cos(/3 — a). In Appendix E we show the branching ratios for all Higgs 
bosons for specific benchmark scenarios to be discussed later. 

7.2 Analysis of Production and Decay Channels 

In the following we study the impact of searches for heavy higgs bosons at ATLAS and CMS. 
To this end, we compute the production cross section and various decay rates for the heavy 
Higgs bosons. We generate the gluon-fusion production cross section at next-to-leading 
order (NLO) using HIGLU [59], taking into account the contributions of the bottom quark 
loop and use the leading order expressions for the partial decay width with the appropriate 
couplings of our model [60, 61]. When relevant, we also consider the vector-boson fusion 
production cross section, using the values quoted in [62, 63]. For charged Higgs production 
we use the NLO results in [64] . In the following we will assume M = Ma = Mh = Mfj +, 
if not specified otherwise, and we discuss in detail the effects of a splitting between the 
neutral and charged Higgs boson masses. 

One of the most interesting channels for the discovery of the pseudoscalar Higgs boson, 
involves the A —)• hZ decay, because the corresponding branching ratio becomes dominant 
for sizable values of cos(/3 —a). There are several experimental studies constraining a{gg —)• 
A) X Br(A —7- hZ), with the light higgs further decaying into bottom quarks [65, 66], tau 
leptons [65] , as well as multi-leptons [49] . 
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CMS: M = (500,600) GeV 
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ATLAS: M = (500,600) GeV 



0.2 0.4 0.6 
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Figure 15: In the left panel we show current exclusion bounds for a{gg —)• A) x Br(A — )• 
hZ —7- Ubb) based on the CMS data [66] . In the right panel we show exclusion bounds for 
a{gg —)■ A) x Br(A —>■ hZ) x Br(/i —)• bb) (blue) and (T{gg —>■ A) x Br(A —)• hZ) x Br(/i —)• 
T~^T~) (green) based on ATLAS data [65]. In both plots we assume equal scalar masses, 
M = 500 GeV (dotted) and M = 600 GeV (dashed), and narrow-width approximation. 
The region below and to the right of the curves is excluded. 


The predictions of our model for both a{gg —)• A) x Br(A —?• hZ) and cr(gg —?• A) x 
Br(A —)• hZ —7- i^i~bb) are presented in Figure 14 in the left and right panels, respectively. 
For the decay rate r(/i —?■ bb), NLO corrections are sizable and therefore we include them 
in our analysis by setting 

T{h ^ bb) = 0.57, (7.9) 

where we use T™ = 4.07 MeV [32] and Bi(Z —>■ ££) = 6.729% for £~ = e~, g,~ [32]. In 
the left panel of Figure 14 we show the contours of cr{gg —>• A) x Br(A —)• hZ) in picobarn 
for 8 TeV proton-proton (pp) collisions in the cos(/3 — a) — tan /3 plane for M = 600 GeV. 
The shape of the contours follows naturally from the fact that the branching ratio scales as 
cos(/3 — a)^, while the production cross section depends only on tan jS. This is no different 
than in a generic two Higgs doublet model [5, 49], but it is particularly relevant in our 
model, since it cannot live close to the decoupling limit, as discussed in Section 4. The 
experimental exclusion bounds from [49] constrain a{gg —)• A) x Br(A —)• hZ) considering 
a multi-lepton final state, but the study is only performed for pseudoscalars with masses 
up to Ma < 360 GeV. 

In the right panel of Figure 14 we show the contours of (T{gg —>■ A) x Br(A —)• hZ —)• 
£~^£~bb) in femtobarn for 8 TeV pp collisions in the cos(/3 — a) — tan/3 plane and M = 600 
GeV. Two branches with suppressed values for a{gg —>■ A) x Br(A —)• hZ —)> ££hb) appear. 
The first branch is the decoupling or alignment limit, where gAhZ vanishes. The second 
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Ma = Mh= 600 GeV, Mh+ = (400,600) GeV 
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Figure 16: In the left panel, we show the exclusion contour for a{gg —>• ^4) xBr(A —)• hZ) x 
Br(/i —)• bb) for M = 600 GeV in the narrow width (NW) approximation (dashed blue) and 
taking into account finite width (FW) effects (solid blue). The black contour additionally 
shows mass splitting effects, assuming Ma = Mh = 600 GeV and = 400 GeV. The 
shaded region inside each contour depicts the excluded area. The right panel shows the 
rescaling factor due to FW effects with respect to the NW approximation extrapolated 
from the CMS analysis [66], for Ma = 500 (600) GeV in pink (green). 


branch is given by the region for which the coupling of the light Higgs h to bottom quarks 
becomes small. 

We consider the measurement of cr{gg —?■ H.) x Br(H —)• hZ) x Br(/i —?■ bh) by ATLAS 
[65] and the measurement of cr{gg —)■ A) x Br(A —)■ hZ —)■ t^l~bb) by CMS [66] with 
= e~ , . Both experiments give their bounds assuming narrow width approximation 

for the heavy scalar. In Figure 15 we compare these bounds (blue curves) from both for 
equal masses of the heavy scalars with M = 500 GeV (dotted) and M = 600 GeV (dashed). 
For both mass choices the ATLAS measurement gives a stronger bound. For M = 500 GeV, 
substantial regions of the model parameter space are ruled out, however for M = 600 GeV 
the model is considerably less constrained. The right panel of Figure 15 also shows the 
ATLAS bounds [65] of a{gg —)• A) x Br(A —)• hZ) with the light Higgs h decaying further 
to tau leptons. The corresponding bounds are shown as green contours for M = 500 GeV 
(dotted) and M = 600 GeV (dashed). These constraints are substantially weaker than the 
corresponding bounds for the h ^ bb decay. 

In the following we consider the impact of finite width effects on the previous bounds. 
In the right panel of Figure 16, we show the rescaling factor for the cross section times 
branching ratio due to finite width effects, extrapolated frm the CMS analysis [66], for 
Ma = 500(600) GeV in pink (green). In the left panel of Figure 16 we first show for 
comparison the exclusion bound from ATLAS data for M = 600 GeV in the narrow width 
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Figure 17: Model predictions for the contours of (T{gg —?■ H) x Br(i7 
for 8 TeV pp collisions and heavy scalar masses M = 600 GeV. 


hh) in picobarn 


approximation. Under the assumption that the scaling effects for ATLAS and CMS are 
similar and assuming sensitivity up to a total width of ~ 100 GeV, we consider finite 
width effects for each point in the cos(/3—a) —tan j3 plane and reinterpret the ATLAS results 
(solid, blue line in the left panel of Figure 16). Although finite width effects significantly 
weaken the exclusion bound, this channel remains the most promising discovery channel 
at the LHC run 11. The bound is further relaxed in the case of a mass splitting, Ma 
M fj+{MH), such that the decay channels A —)• H^W~{HZ) open up. Our discussion in 
Section 6 showed that such a mass splitting is only allowed between the pseudoscalar and 
the charged Higgs boson. We present the bound for a{gg —>■ A) x Br(A —)• hZ) x Br(/i —)• bb) 
including finite width effects for Ma = Mh = 600 GeV and Mff+ = 400 GeV in the left 
panel of Figure 16 (black). 

In the following we will consider the experimental bounds from searches for the neutral 
CP-even Higgs boson H. There are two channels of particular interest, the CP even scalar 
decaying into light Higgs bosons H ^ hh and the CP even scalar decaying to vector bosons 
H ^VV with V = W,Z. 

In Figure 17 we present predictions for a{^gg —)■ H) x Br(H —)■ hh) in picobarn for 
8 TeV pp collisions in the cos(/I — a) — tan /3 plane for M = 600 GeV. From (7.4) we 
observe that the self coupling gnhh is proportional to cos(/3 — a) and has an explicit Ma 
dependence. For cos(/3 — a) > 0 we observe two branches of contours with suppressed 
a X Br. The first branch approaches zero at cos(;0 — a) = 0, and for the second branch 
both the coupling gnhh and the production cross section become small. Predictions for 
(^{99 H) X Br(iL —7- hh) are comparable to the ones in a generic two Higgs doublet 
model of type H [5]. Similar to the pseudoscalar case, the experimental exclusion bounds 
for a{gg — H) x Br(i7 — hh) [49] are only available up to Mh < 360 GeV. However for 
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M = 600 GeV M = 600 GeV 




Cos(j6-a) Cos(j6-a) 

Figure 18: Contours of a{pp H + X) x Br{H —)■ VV)/{a{pp H + X) x Br{H —)• 
^^))sM (right panel) and a{pp —> qqH) xBr(i/ —)■ VV)/{cr{pp —> qqH) ^ FF))sm 

(left panel). The heavy scalar masses are set to M = 600 GeV. 


the CP even Higgs, the model predictions seem to be much below the present experimental 
sensitivity. 

The most important search channel for the heavy CP even neutral Higgs boson H 
is the inclusive production with subsequent decay oi H ^ VV with V = W, Z. In our 
specific model, in particular, there is an interesting region of parameter space in which 
the vector boson fusion production is competitive with the gluon fusion production due 
to the behavior of k^. Normalized to the corresponding SM Higgs production and decay 
processes for a SM Higgs of mass Mjj, we have for gluon fusion and vector boson fusion 
production processes, respectively. 


a(gg ff) x Br(ff VV) 

(a(gg ^ H) x Bt{H ^ VV))sm 
a{pp qqH) x Br{H VV) 

{a{pp qqH) x Br{H W))sM 



(7.10) 

(7.11) 


where denotes the correction from a bottom quark in gluon fusion with respect to the 
leading top contribution. We take the SM total width for a heavy Higgs of mass Mh 
from the LHC Higgs Cross Section Working Group [62, 63, 67]. 

In Figure 18 we present theoretical predictions for contours of inclusive heavy neutral 
CP even Higgs production (left panel) and vector boson fusion production (right panel) with 
subsequent decay into H —> VV, using (7.10) and (7.11), for M = Ma = Mh = Mh+ = 
600 GeV. The vector boson fusion is governed by Ky and becomes strongly suppressed for 
small cos(/3 —a). The gluon fusion production mode in (7.10) is suppressed for small values 
of or for small Ky and this effect shows in the inclusive production mode above. We 
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M = (500,600) GeV 



Cos(j6-a) 

Figure 19: Exclusion bounds for a{pp H + X) x Bic{H — )■ VV)/{a{pp —?■ H + X) x 
Br(F ^ FF))sm of CMS [68] for M = 500 GeV (dotted) and M = 600 GeV (dashed). 

observe that for small both production cross sections become competitive. The theory 
prediction for these two observables differs from a two Higgs doublet model of type II only 
by the different scaling of the width Tjif and the contribution of the bottom quark to gluon 
fusion, which is small for tan /3 rsj 0(1). 

The CMS collaboration has reported updated results from an inclusive search for a 
heavy Higgs decaying into W~^W~ and ZZ in the range of Mh = 145 — 1000 GeV [68]. 
They consider both fully leptonic and semileptonic final states. In Figure 19 we illustrate 
those bounds for M = Ma = Mh = Mh+ with M = 500 GeV (dotted) and M = 600 GeV 
(solid). We observe that this search mode is competitive with the bounds obtained from 
the A —>■ hZ channel. We note that for the neutral CP even Higgs analysis no finite width 
effects have been taken into account, although we expect sizable finite width effects in a 
large region of parameter space, compare the left panel of Figure 13 . 

The CMS collaboration also performed an analysis for a heavy neutral Higgs boson 
decaying into W~^W~ in vector boson fusion production channel in the mass range Mh = 
110 — 600 GeV [69]. The observed signal significance is close to the SM prediction for a 
Higgs of Mh = 300 — 600 GeV, and hence from the right panel of Figure 18 it follows that 
there is no sensitivity to the preferred parameter region from this search. 

Searches for heavy charged Higgses have been performed by both ATLAS and CMS 
collaborations. In particular, they searched for production modes in association with a 
single top, a{bg —)■ H~t), or top and bottom quarks, cr{gg —)■ H~tb), with subsequent 
decays into third generation fermions: H~ —)• tb and H~ —)• Ti'r [70-72]. The most recent 
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limits are 

—)• Ti/) < 0.153 pb — 0.026 pb for Mj:^+= 300 — 600 GeV , (7.12) 

Bi{H~ ^ tb) < 6pb —4pb for Mjj+ = 300 — 600 GeV, (7.13) 

assuming J3r(H~ —)• tv) = 100% and —)• tb) = 100%, respectively. These val¬ 

ues are below the expected production cross section, a{pp —)■ H~t{b)) = 70fb — 6fb for 
Mfj+ = 300 GeV — 600 GeV and tan /3 « 2 (lower values of the production cross section 
occur for 2 < tan/3 < 6) [64]. A heavy charged Higgs boson is therefore not constrained 
for the parameter region of interest, through current direct search limits. 

For a heavy charged Higgs Mh+ k, 360 — 400 GeV,cos(/3 — a) > 0.3(0.2) and tan/3 = 
2(4), the decay channel H'^ —)• hW~^ dominates over —>■ tb. The branching ratio can 

become as large as Br(F7’^ —>■ hW^) « 85% for tan/3 = 2.5, cos(/3 — a) = 0.6. For a lighter 
charged Higgs, this is slightly less pronounced and we find Br(i7'’' —)• hW~^) ~ 70% for 
tan/3 = 2.5, cos(/3 — a) = 0.6 and Mfj+ = 400 GeV. 


8 Origin of the Effective Yukawa Couplings 

In this section we present an example of the origin of the effective Yukawa couplings at the 
TeV scale for the bottom quark sector. Similar considerations can explain the generation 
of the other effective light quark Yukawa couplings in our model. A complete description 
of the UV completion is beyond the scope of this paper. 

A possible completion of the Froggatt Nielsen model may introduce new colored vector¬ 
like fermions or additional scalar doublets [73], whose masses determine the suppression 
scale A in the expansion parameter (2.5). Since in our model the flavor breaking scale is 
identified with the electroweak scale and the expansion parameter is fixed by the ratio of 
bottom and top quark masses e = mb/mt, the UV scale is constrained to be of the order 
of A ~ 1 TeV. 

The relevant operators that would provide a UV completion for the bottom Yukawa 
interactions are 


Tluv = yi bL Hd r]R y2 t]r Hu 'tpL + Vs tpL Hd bn + MnruVR + , (8.1) 

such that after integration of the heavy fields the effective Lagrangian is given by 

C-EFT = bi Hd bR , ( 8 . 2 ) 

with 




yi y 2 y3 VuVd 

MrjM^ 2 


(8.3) 


The corresponding diagram is given in Figure 20 in which the new vector-like fermions 
carry quantum numbers 


yL,R ~ (3, 1, -1/3, 2), ^|;L,R ~ (3, 2, 1/6,1) 


(8.4) 
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Figure 20: Diagram in the full theory, which generates the Yukawa coupling between the 
Higgs and the bottom quarks after integrating out the heavy vector-like fermions 'll:, t]. 


with respect to the groups {SU{3)c, SU{2 )l, U{1)y, f7(l)i?). 

From (8.3) is follows that for fixed yi = 2/2 = 2/3 = 1 and y'^ G [0.5,1.5] this predicts 
the masses = M,p ~ A = 1 TeV. It is evident that slightly larger fundamental Yukawa 
couplings yi, y 2 and 2/3, allow for heavier vector-like fermions, while any tan/3 /:?> 1 or 
tan/3 <C 1 lead to lower mass scales. In the spirit of avoiding hierarchies between the 
fundamental couplings, including the top Yukawa coupling, we shall consider the ratio 
Vi/yt 0(1) with i = 1, 2, 3. This constrains the masses of the vector-like fermions to be 
at most of the order of a few TeV. In particular, we define a generic mass M = 
and an average fundamental Yukawa coupling y = ( 2 / 12 / 22 / 3 )^^^) such that 


= 


ty 


y3 y2 

y<^ 2s 1 + 3 ^ 


(8.5) 


In Figure 21, we show the expected masses of the new fermions for varying tan/3 and 



Figure 21: Masses of the new fermions in the UV completion depending on tan /3 and for 
three different values of the average Yukawa coupling y = 1,1.5, 2 (from bottom to top). 
Fermion masses below the solid red line are excluded by current LHC bounds, while the 
dashed red line shows the expected exclusion reach for the 14 TeV run of the LHC. 
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fixed = 1, for three different values of average Yukawa couplings y = 1,1.5, 2 (from 
bottom to top). These predictions for the expected masses remain the same for y = 1 and 
change at most by 15% (25%) for y = 1.5 (2) for the hrst generation quarks and at most 
10% (20%) for second generation quarks. 

The solid and dashed red lines in Figure 21 indicate the present and projected experi¬ 
mental bounds from searches for pair produced heavy quarks at the LHC. These searches 
have been performed both by ATLAS and CMS, and exclude vector resonances with masses 
of 600 — 800 GeV [74-76], depending on the decay mode, with some channels already prob¬ 
ing top partners T up to 900 GeV for Br(r —)• W~^b) = 100% [77]. The next run of the LHC 
has a projected reach of M > 1.2 (1.4) TeV for 20 fb“^(100 fb“^) and Br(r —>■ W~^b) = 50% 
[78]. Searches for heavy vector-like quarks in single production have also been considered 
[79-81] and could be much more effective as a discovery channel for sufficiently heavy 
vector-like quarks compared to the previously mentioned pair production searches. How¬ 
ever, the LHC reach in the single production channel depends very strongly on the model 
parameters which define the couplings of the heavy quarks to SM quarks. A reinterpreta¬ 
tion of any of the existing LHC bounds in single heavy quark production channels would 
demand a detailed study of production cross sections and decay branching ratios for a spe¬ 
cific UV completion. Similarly, a specific UV completion would be subject to constraints 
from electroweak precision measurements as well as from flavor physics [73]. The latter 
have been addressed in some detail in the original Giudice-Lebedev paper [3]. 

9 Benchmark Scenarios 

The global fit to Higgs signal strength measurements discussed in Section 3 universally 
constrains the allowed parameter space to two branches within cos(/3 — a) = 0.35 — 0.8. 
Smaller values of cos(/3 — a) < 0.35 are in principle possible for tan/3 > 5, but such large 
values of tan/3 are in tension with flavor observables. Electroweak precision observables 
and collider searches for the extra scalars provide additional constraints that narrow the 
parameter space significantly. In the following we examine the allowed window in the 
cos(/3 — a) — tan/3 plane and specify three benchmark points, that highlight the interest¬ 
ing features for the phenomenology of this model, and for which we give a detailed list of 
couplings, production cross sections and decay widths. 

As a result of the discussion in Sections 6, the combination of constraints from fla¬ 
vor physics, electroweak precision observables, unitarity and perturbativity lead to a con¬ 
strained region of mass values for the additional Higgs bosons ~ Mh ~ 500 — 600 
GeV and Mh+ ~ 360 — 500 GeV. Perturbativity puts an upper bound of 600 GeV on 
the neutral Higgs masses and requires a splitting between the neutral and charged Higgs 
masses of Ma,h — Mh+ ^ 100 GeV. Electroweak precision measurements exclude the left 
branch of the global fit to Higgs coupling measurements for values of tan /3 < 4.5. In addi¬ 
tion, tree-level contributions to meson-antimeson mixing put an upper bound of tan /3 < 5, 
while loop contributions from charged Higgs exchange result in a lower bound tan /3 > 1.5. 
Collider searches for the two heavy neutral Higgs bosons further constrain the allowed pa- 
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rameter space and probe the right branch of the global Higgs fit for cos(/3 — a) = 0(0.5) 
and tan/3 <3. As a result, there is a specific window of allowed masses as well as values of 
cos(/3 — a) and tan /3, which translates into a precise prediction for searches for the extra 
scalars and constrain the possible deviations in the SM Higgs couplings. In Figure 22, we il¬ 
lustrate this window by showing the 95% CL region of the global ht to ATLAS Higgs signal 
strengths measurements (red shaded area), the region preferred by electroweak precision 
constraints (shaded green) and the bound induced from flavor constraints (solid purple 
contour), as shown in Figure 9. Further, we superimpose the bounds derived from the 
ATLAS and CMS measurements of (T{gg —)• A) x Hr (A —)■ hZ —)• (solid blue) and 

cr{pp H + X) X Br{H —)• VV) (solid orange). In the left panel we assume scalar masses 
of Mh = Ma = 600 GeV, Mjj+ = 450 GeV, and in the right panel Mh = Ma = 500 GeV, 
Mfj+ = 360 GeV. The gray shaded area is excluded, the overlap of the light green and red 
regions is allowed. 

Comparing the two plots in Figure 22, bounds from flavor physics as well as collider 
constraints become weaker for larger masses. The area in agreement with electroweak 
precision bounds is slightly larger for smaller mass splittings, but similar for the two ex¬ 
amples given in Figure 22. The right boundary of the right branch of the global Higgs 
fit is close to the contour of Kb = —1, for which the Higgs coupling to bottom quarks has 
the same size, but opposite sign compared to the SM one. The left boundary of the right 
branch is close to Kb = —0.5. For all of the allowed parameter space, we can therefore infer 
“1 ^ ^ —0.5. In addition to the sign and the reduction of the Higgs bottom coupling, 

we find a universal enhancement of the Higgs charm couplings. Both can in principle be 
probed by measurements of exclusive radiative Higgs boson decays, which can test the sign 
of Kb at the 14 TeV LHC, and establish possible departures from the SM Higgs charm 
couplings of the order of 20% at a prospective 100 TeV collider [82, 83]. In the presence 
of a Higgs portal to dark matter, such corrections to the Higgs couplings to quarks could 
signihcantly modify the direct detection cross section [84]. 

In Table 4 and 5 we give the values for the Higgs couplings, signal strengths, production 
cross sections and branching ratios for three representative benchmark points indicated by 
black crosses in Figure 22. Typical values of cos(/3 — a) ss 0.4 — 0.55 and tan /3 ss 3 — 4.5 
are considered. In all cases, Kt ^ 1, implying a gluon fusion production rate of order of the 
SM one. 

Benchmarks la and lb allow for larger values Ma,h ~ 600 GeV and a charged Higgs 
mass Mfj+ ~ 450 GeV, close to the 2cr bound derived from the experimental b ^ s'j 
measurement in a type H two Higgs doublet model with tan /3 > 2. 

In Benchmark la, the tree-level gauge boson and down type fermion third generation 
couplings are suppressed by factors of order 20% and 40%, respectively, while the Higgs 
coupling to charm is enhanced by about 20%. The sizable suppression of Kb yields a 
suppression of the branching ratio into gauge bosons and hence of the corresponding signal 
strength of those channels. The charm signal strength instead, is increased by a factor 
~ 2 — 3 (depending on the production mode) due to the combined effects of an enhancement 
in Kc and a suppression in Kb and Ky- All other vector boson fusion and VH production 
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Benchmark 1 : Ma = Mh = 600 GeV, Mfj+ = 450 GeV, 
la cos(/3 — a) = 0.55 , tan P = 3, 
lb cos{(3 — a) = 0.42 , tan/3 = 4.5, 

Light Higgs Gouplings; 


la Kt = 1.02 , Kv = 0.84, Kh = Kr = —0.61, Kc = 1.22 , Ks = —0.41 
lb Kt = 1.00 , Kv = 0.91 , Kh = Kr = —0.96 , Kc = 1.02 , Ks = —0.95 

Higgs Signal Strength: 


la 

flV 

P-y 

Pb 

l^c 

lb 

PV 

Py 

Pb 

Pc 

^gg^h 

1.38 

1.21 

0.74 

2.95 

^gg^h 

0.96 

0.91 

1.09 

1.22 


1.33 

1.17 

0.71 

2.84 

^tt^h 

0.90 

0.85 

1.02 

1.14 

O'VBF, O'VH 

0.89 

0.78 

0.48 

1.91 

0'VBF,(yvH 

0.74 

0.70 

0.84 

0.94 


Heavy Scalar Production Cross Sections for la (lb): 


8 TeV: cr{gg A) = 78(36) fb, <y{gg —> H) = 32(21) fb, 
a{pp —^ H~t{b)) = 9(4) fb , 

14 TeV: a{gg ^ A) = 361(157) fb, a{gg ^ H) = 166(97) fb, 
(t{pp —)■ H~t{b)) = 63(25) fb , 


Heavy Scalar Decay Modes: 



Total Width for la (lb): 

Pfe = 2.22(3.71) MeV, Pa = 24.6 (16.3) GeV, P^ = 36.4 (26.1) GeV , 

Th+ = 10.2(5.8) GeV. 

Table 4: Values for the Higgs signal strength, heavy scalar production cross sections for 
the dominant channels at the LHC, partial and total widths for the benchmarks la and 

lb. 


channels are suppressed with respect to the SM, in particular the /i —)■ 66 search mode. 

In Benchmark lb all tree-level fermion and gauge Higgs couplings are within less than 
5 — 10% of the SM expectations, hence the signal strengths in gluon fusion production are 
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Benchmark 2 : Ma = Mh = 500 GeV, Mj^+ = 360 GeV, 
cos(/3 — a) = 0.45 , tan /3 = 4, 

Light Higgs Gouplings: 

lb Kt = 1.01 , Ky = 0.9 , Kb = Kr = —0.81 , Kc = 1.1 , Kg = —0.71 
Higgs Signal Strength; 


2 

Pv 

P-y 

9b 


^gg^h 

1.15 

1.07 

0.94 

1.76 


1.09 

1.02 

0.90 

1.67 

O'VBF, O'VH 

0.86 

0.80 

0.71 

1.32 


Heavy Scalar Production Cross Sections; 


8 TeV; a{gg A) = 130 fb, 
14 TeV; a{gg A) = 546 fb, 

Heavy Scalar Decay Modes; 


<^{99 —> 77) = 53 fb , a{pp —)■ H t{b)) = 12 fb , 
(T{gg H) = 224 fb, a{pp —^ H~t{b)) = 66 fb, 


A 

Fi/TA 

Zh 

56.6% 

W-H+ 

23.3% 

bb 

5.3% 

tt 

12.4% 

T+T~ 

0.66% 

tc 

1.1% 


H 

F^/FH 



WW 

45.4% 

77+ 

r^/FH^ 

zz 

21.8% 

hW 

71.8% 

hh 

11.5% 

tb 

27.8% 

W-H+ 

12.6% 

TU 

0.4% 

tt 

3.65% 




Total Width; 

P/j = 3 MeV, Pa = 10.7 GeV, Fh = 15.7 GeV, r//+ = 3 GeV. 


Table 5; Values for the Higgs signal strength, heavy scalar production cross sections for 
the dominant channels at the LHC, partial and total widths for the benchmark 2. 


also within 5 — 10% of the SM ones, with the exception of a 20% enhancement in pc- All 
vector boson fusion/VH production channels are suppressed with a maximal suppression 
of 25 — 30% in the case of light Higgs decaying into gauge bosons. 

Benchmark 2 allows for the smallest possible values of Mjj+ = 360 GeV compatible 
with the 3a bounds derived from the experimental b ^ measurement in a type H 2HDM 
with tan/? > 2. Benchmark 2 has a similar tendency in the couplings of gauge bosons and 
fermions to the light Higgs boson and hence in the corresponding signal strengths as in 
Benchmark la, but with percentual effects in the deviations from SM predictions that are 
a factor 2 — 3 smaller. In addition to improving signal strength measurements, the ongoing 
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Ma = Mh = 600 GeV, M^+ = 450 GeV 
6 


0.4 0.6 

Cos(yS-Q') 


0.4 0.6 

Cos(jS-a) 


Ma = Mh = 500 GeV, Mh+ = 360 GeV 
6 


Figure 22: Summary plots showing constraints from flavor observables (purple con¬ 
tour) and direct collider searches for A —>■ hZ —)• (blue contour) as well as 

H —7> W^W~jZZ (orange contour), where the gray shaded area shows exclusion. The 
red shaded region is allowed at the 95% CL from the global fit to Higgs signal strength 
measurements to ATLAS data. The green area highlights the allowed region from elec- 
troweak precision observables, perturbativity and unitarity constraints. The panels corre¬ 
spond to Ma = Mh = 600 GeV and Mfj+ = 450 GeV (left), and Ma = Mh = 500 GeV 
and Mfj+ = 360 GeV (right). The black crosses in both panels indicate the benchmark 
scenarios. 


run of the LHG will probe these benchmarks by direct searches for the additional Higgs 
bosons. All three benchmark scenarios will be primarily tested by the search for A —)• Zh 
and H —)• VV, that have branching ratios of 55% — 75%, depending on the scenario. In 
the case of ff —)■ VV, the inclusive and vector boson fusion production modes will play 
a complementary, relevant role. In addition to these discovery channels, other interesting 
search modes such as A, H —)• W^H~, H —)■ hh, A —)• tt, —>■ hW~^, and 77+ —)> tb 

would yield additional valuable information about this model. The mass splitting between 
neutral and charged scalars give rise to an additional decay chain, that can potentially 
allow to discover the charged Higgs even for masses of Mfj+ ss 360 — 400 GeV, in particular 
for the subsequent decay of 77+ —)• W~^h. Although challenging due to the small branching 
ratio, a novel channel in these scenarios is A —)• tc. 

Predictions for particular observables can be computed from the information provided 
in Table 4 and Table 5. Finite width effects play a relevant role and in the case of A —?■ hZ 
we have compiled them in the right panel in Figure 16. 

Finally, improved measurements of flavor observables, in particular in the neutral 
system could additionally constrain the parameter space significantly. 
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10 Conclusion 


In this article we propose an explanation to the hierarchies in fermion masses and mixings 
based on a Froggatt-Nielsen mechanism, in which two Higgs doublets play the role of the 
flavon. Therefore, the underlying flavor symmetry is broken at the electroweak scale. The 
flavor charges are fixed to reproduce the SM quark mass hierarchies and CKM mixing 
angles up to rescalings, that have no effect on any physical quantity. As a result, this two 
Higgs doublet flavor model can be described by few effective parameters, the masses of the 
extra scalars Mh, Ma, Mu+, cos(/3 — a) and tan/3. This allows us to present our main 
findings in the cos(/3 — a) — tan/3 plane for fixed mass values, as shown in Figure 22. 

Modified interactions between the SM-like Higgs h and quarks are characteristic for our 
two Higgs doublet flavor model, leading to strong constraints from Higgs signal strength 
measurements. The results of our Higgs global fit to ATLAS and CMS data constrain 
possible deviations of the couplings of the light Higgs to fermions and gauge bosons with 
respect to the SM ones, and select sizable values of cos(/3 — a) 0(0.5). This implies a 
suppression of the tree-level couplings of the Higgs to gauge bosons, which is proportional to 
sin(/3 —a) as in any two Higgs doublet model and therefore a suppressed vector boson fusion 
production rate with respect to the SM. The alignment/decoupling limit cos(/3 — a) = 0 is 
excluded for all values of tan /3, since in this limit our model approaches the Babu-Nandi- 
Giudice-Lebedev model for which there is a factor of three enhancement for coupling of 
the light Higgs to bottom quarks. The Higgs global fit allows for two branches in the 
cos(/3 — a) — tan/3 plane (red shaded areas in Figure 22) with opposite sign of the bottom 
Yukawa coupling. However, other constraints end up singling out the branch with values 
of the SM normalized light Higgs-bottom Yukawa coupling between —0.5 and —1. On this 
branch the light Higgs-top Yukawa coupling is close to its SM value, implying gluon fusion 
signal strengths of 0(1). Furthermore, on this branch, the coupling of the light Higgs to 
charm quarks is universally enhanced by up to 30%, leading to a possible enhancement 
of the Higgs to charm signal strength by a factor of three. Both the negative sign of 
the bottom Higgs coupling as well as the enhanced Higgs to charm signal strength can in 
principle be measured at a high luminosity/energy collider through exclusive Higgs decays 
with a final state photon, such as /i —?• Ty and h —?■ J/ifj'y. 

Flavor changing neutral currents arise at tree-level, mediated by the light Higgs as well 
as the extra neutral scalars. Remarkably, light Higgs FCNCs become automatically small 
for the branch of the global Higgs fit with negative light Higgs-bottom Yukawa coupling. 
While the masses of the extra neutral scalars are constrained to be larger than 500 GeV in 
this region, we need a mild fine-tuning of 0(10%) in the Yukawa couplings in order not to 
exceed the strongest constraint from Bd — Bd mixing (shown as purple contour in Figure 22). 
These tree-level FGNGs result in an upper bound of tan/3 <5.5. Moreover, contributions 
from box diagrams with charged Higgs exchange can compete with the tree-level diagrams 
for low tan/3 and exclude values of tan/3 < 1. Thus the interplay of tree-level and loop 
contributions in flavor observables predicts 5.5 > tan/3 > 1. Interestingly, if we discard the 
explanation of the GKM angles by the two Higgs doublet flavor model, we find almost no 
constraints from flavor observables in the region preferred by the global Higgs fit. As in any 


- 42 - 


two Higgs doublet model, charged Higgs exchanges also induce FCNCs through penguin 
diagrams, for example b ^ s'j, which imply a lower bound on the charged Higgs mass of 
360 GeV for tan/3 > 2. 

The two Higgs doublet flavor model offers exciting possibilities for direct collider 
searches for the additional Higgs bosons. Electroweak precision observables, perturba- 
tivity and unitarity constraints choose a preferred range of masses and mass splittings for 
the new heavy scalars. In particular, almost degenerate values for the CP-odd and CP-even 
Higgs boson masses and sizable splitting between the neutral and charged Higgs masses 
are strongly favoured. This opens the opportunity of new decay channels, A —)• H^W~ 
and H —)■ H^W~ in addition to the regular decay channels H —)• W~^W~jZZ, A —>■ hZ, 
that are importantly enhanced in the cos(/3 — a) ~ 0(0.5) region. The latter are the 
leading discovery modes for these scalars (present bounds are shown by blue and orange 
contours in Figure 22). Furthermore, the cos(/3 — a) dependence of the H W^W~, HZZ 
couplings are of particular relevance because the vector boson fusion production mode can 
compensate for the suppression of the gluon fusion production mode of the CP even Higgs 
in the relevant regions of parameter space. Direct searches for a charged Higgs boson are 
not sensitive for masses compatible with the flavor constraints, however future searches via 
Higgs decay chains with the subsequent decay —)• W~^h may be promising. The other 

possible decay of heavy Higgs bosons to the SM Higgs is in the channel H ^ hh with 
branching ratios of order 10%. 

The fact that the flavor symmetry is broken at the electroweak scale predicts a UV 
completion in the few TeV range, as well as a low value of tan /3 in agreement with flavor 
constraints. The necessity of new physics at the TeV scale provides an additional motivation 
for the search for new vector-like fermions at the run H of LHC. 

We conclude, that in the two Higgs flavor model constraints from flavor observables, 
Higgs precision measurements, direct heavy Higgs searches, and precision electroweak ob¬ 
servables, as well as unitarity and perturbativity restrictions on the theory, can be fulfilled 
simultaneously. We propose three benchmark scenarios in this region, that highlight dif¬ 
ferent characteristics of the two Higgs doublet flavor model (black crosses in Figure 22). In 
Table 4 and 5 we provide all the relevant information to compute production cross sections 
and decay rates for these benchmark scenarios and test the two Higgs doublet flavor model 
at the run H of LHC. 
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A The Higgs Potential 


In this appendix we consider the scalar potential and related topics. 

The fact that HuHd carries a flavor charge strongly constrains the scalar potential. 
We need a (soft) source of flavor breaking in order to generate a 6 -term. We consider this 
additional source of flavor breaking to be irrelevant for the texture of the Yukawa couplings. 
The potential reads then 


V{Hu, Hd) = 111 Hi Hu + Hi Hd - [6 Hu Hd + h.c] (A.l) 

+ ^ {Hi Hu? + ^{hI Hd? + Ag {Hi Hu){hI Hd) + ^{Hi Hd){Hl Hu) , 

in which HuHd = H?{ia 2 )Hd- Note that the potential is the same as in a generic CP- 
conserving two Higgs doublet model, see for example [9, 12], whith As = Ae = A 7 = 0. 

In order to diagonalize the potential, we introduce the neutral scalar mass eigenstates, 


(h\_(ca -Sa\ /Rei7°\ 

\h) ~ \sa C„ ) [ReHUJ ’ 


(A.2) 


with the mixing angles Cq = cos a and Sa = sin a as well as the pseudo-scalar mass 
eigenstates 


/ S/s -c^\ /lm7l0\ 
Sfs ) [imH^J ^ 

and similarly for the charged mass eigenstates H^, 

[h-J ■ 


(A.3) 


(A.4) 


Performing these rotations the explicit formulas for masses of scalar fields can be obtained, 
see for details for example [9, 12]. 

Einally from the scalar potential we obtain all couplings between the scalars [5, 9, 12]. 
In particular, relevant for our analysis are the coupling between the heavy scalar H and the 
light Higgs h given in equation (7.4), and the coupling of the light Higgs h to two charged 
Higgses H^ 


9hH+H- 


{2Ma - 2M^± - + 2(M^ - m^,) 


2 n 


■S 2/3 


(A.5) 
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B Box Diagrams and Loop Functions 


In this appendix, we collect the contributions to the Wilson coefficients (5.1) from box 
diagrams and the relevant loop functions [37, 42, 85]. For K — K mixing we have the 
following Wilson coefficients: 


C 


sd 

1 ,box 


ml 


w 


Cl ,box + C^,hox + Kd Cl^box) 


Ibvr^ 

with At = VtdVtl, Ac = Fed Fes and 

C'gbox = + xl) D2{mt, mw) - 8 x? Do{mt, mw) 

2x^ 

+ ~7f [txiw D2{mt,mw,MH±) - 4m^ Do(mt, mw, 
i/3 

3,2 

+ TT C>2{nxt, Mu±) , 

i/3 

C'gbox = (4xc + Xe) mw D2{mc, mw) - 8 x^mw Do{mc, mw ), 

C'gbox = (4xct + x^i) m^ D2{mc, mt, mw) - 8 x^j m^ Do{mc, mt, mw) 
2x^ 


x2 

+ '^w D2{mcmt, Mh±) 


(B.l) 


(B.2) 

(B.3) 


+ [m\^D2{mc,mt,mw,MH±) - Am'^ Do{mc,mt,mw, Mh±)] 

i/3 


(B.4) 


in 


which xt = m^jm^^Xc = m^/mw and Xct = mcmt/m‘^. For mixing, we 


have 


C 


bq _ G'iml 


:oq _ ^ t •••'W ( \t \2 (~it 
1 ,box 167 r 2 ^'^bq) '-^l,box’ 


c: 


bq _ 

2 ,box 


Gpmyy 4mg , 2 ^ 

,2 y'^bq) *-^2,box, 


IGvr^ 


m 


W 


c: 


bq 
4 ,box 


_ Gj.m'^ 2mqmb , ^ .2 


IGvr^ 


m 


w 


i'^bq) ,box , 


with At = V^^ytq and {q = s, d) and 


G- 


1 ,box 


c: 


bq _ 


2 ,box 


G'p mw 

2 2 
rriqUif^ 

IGvr^ 


Gp 

Ami /■ 

IGtt^ 

™2 F 
i'V 

G\ m^ 

rriqmb 

IGvr^ 



(A* 6 ,)"Ci,box, (B.5) 


(A*6,)2C2,box, (B.G) 


^5 ,box ? 


(B.7) 


F*! ,box IXlw 


mt 


^/3T7^^2(mt,Mt:t±) + i/3 D2{mt,mw,MH±) 


(^ 2 ,box = Xt mw [F>o(?nt, Aftd±) + 2Do{mt,mw,MH±)] 


Q.box = Xt mw 4iAo(mt,Mt:t±) + [tj + ^] D2{mt,mw, Mh±) 
L ^ i/3 ^ 

- 4t^ Xt mwD2{mt,mw,MH±) , 

Cb.box = Xt mw [F>2(nrt,Mtit±) + 2F)2(m3, w-w, . 


(B. 8 ) 

(B.9) 


(B.IO) 

(B.ll) 
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The loop functions are given by 


Do(mi,m2,Mi,M2) = 


m‘ 




(m| — — M|) 


+ 


+ 




{Ml — mf )(Ml — 1712 ){Ml - Ml) 


A^|log(S) 




D2{mi,m2,Mi,M2) = 


m< 




(m| — m\){ml — Ml){m\ — M|) 


+ 


+ 




{Ml — m\ ){Ml — 1712 ){Ml - Ml) 

M|log(^) 


and for i = 1, 2 


Di{mi, Ml, M 2 ) = lim Di{mi,m 2 , Mi, M 2 ), 

m2^mi 

Di{mi,Mi)= lim Di{mi, Mi, M 2 ), 

M2^Mi 

D 2 {mi, Ml, M2) = D 2 {mi,Mi,M 2 ) - £> 2 ( 0 , Mi, M2 ). 


(B.12) 


(B.13) 


(B.14) 

(B.15) 

(B.16) 


C Random Parameter Generation and Running 

In order to find sample parameter points, we generate random fundamental Yukawa cou¬ 
plings with and \ylf\ G [0.5,1.5] and G [0, 27r]. The effective Yukawa 

couplings (2.7) have to reproduce the quark masses and Wolfenstein parameters in Table 
6 in Appendix D. To this end we perform a fit; with symmetrized 2cr errors and require 
X^ < 10. 

In order to obtain the new contributions to K — K and — Bs,d mixing we compute 
the Wilson coefficients with these effective Yukawas, including the tree-level (5.2) and 
one loop Wilson coefficients given by Appendix B. These Wilson coefficients are at the 
high scale /x = ruh and y = Mh,M^±,Ma, respectively. The next step is running the 
Wilson coefficients from the electroweak scale to the scale at which the matrix elements 
are evaluated, /x = 2 GeV in the case oi K — K mixing and y = mf, in the case of Bg^d — Bs,d 
mixing. For K — K mixing we use [41] 

{K\'Hti=‘^\K), = Y,Y. ^''^Cf{y)Bl^{K\Qf\K), (C.l) 

j=l r=l 


- 46 - 








in which rj = asifj) / as{mt), ai,b^p''^ and are ’’magic numbers” collected in [39] and 
are the B parameters collected in Table 8. The matrix elements are given by 

{K\Qf\K) = ^MKfi, 

{K\Qf\K) =Nr( ) MKfh (C.2) 

\md + msj 

with Nr = (—5/24,1/24,1/4,1/12) for r = (2, 3,4, 5) and Mk and md + rus again given in 
Table 8. For Bd^s ~ Bd^s mixing, (C.l) and (C.2) hold with the obvious replacements. The 
corresponding ’’magic numbers” can be found in [41], and all other parameters in Table 9. 

D Numerical Input 

In this Appendix we collect the numerical input used throughout this paper. 


Quark Masses in GeV [86] Wolfenstein Parameters [87] 


mu{mz) 

0.00127 

±0.0005 

A 

0.22551 

±0.00091 

rridimz) 

0.0029 

±0.0012 

A 

0.813 

±0.035 

nisimz) 

0.055 

±0.016 

fj 

0.342 

±0.024 

rucimz) 

0.619 

±0.084 

P 

0.149 

±0.033 

rribimz) 

2.89 

±0.09 




mt{mz) 

171.7 

±3.0 





Table 6: Quark masses and Wolfenstein parameters at the electroweak scale. Errors are 
symmetrized. 


Couplings and Boson Masses [32, 88] 

ae{mz) 

1/127.9 

as{mz) 

0.1185 ±0.0006 

mz 

91.1876 ±0.0021 GeV 

mw 

80.385 ± 0.015 GeV 

Gf 

1.16638-10-5 GeV-2 


Table 7: Gauge boson masses and couplings. 
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Parameters va. K — K mixing [32, 89] 


Bf 

0.537 ±0.007 ±0.024 


0.620 ±0.004 ±0.031 

B^ 

0.433 ±0.003 ±0.019 

B^ 

1.081 ±0.006 ±0.048 

B^ 

0.853 ± 0.006 ± 0.049 

fk 

156.2 ±0.2 ±0.6 ±0.3 MeV 

Mk 

497.614 ± 0.024 MeV 

rUg + rud 

135 ± 18 MeV 


Table 8: Parameters relevant for K — K mixing. 


Parameters in mixing [90-92] 


Bf 

0.85 ±0.03 ±0.02 

Bi 

0.73 ±0.03 ±0.01 

Bi 

0.88 ±0.12 ±0.06 

Bf 

0.95 ±0.04 ±0.03 

Bi 

1.47 ±0.08 ±0.09 

fSd 

186 ± 4 MeV 

Mb, 

5.27942 ± 0.00012 GeV 

nib ± md 

4.29 ± 0.09 ± 0.08 ± 0.02 GeV 


Parameters in Bg — Bg mixing [90-92] 


Bf 

0.86 ±0.03 ±0.01 

BI 

0.73 ±0.03 ±0.01 

BI 

0.89 ±0.10 ±0.07 

BI 

0.93 ±0.04 ±0.01 

BI 

1.57 ±0.07 ±0.08 

IBs 

224 ± 5 MeV 

Mb. 

5.36668 ± 0.00024 GeV 

mb ± mg 

4.38 ± 0.09 ± 0.08 ± 0.02 GeV 


Table 9: Parameters relevant for Bfi g — B^ g mixing. 
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E Branching Ratios 




Table 10: Branching ratios as a function of cos(/3 — a) for the light neutral scalar (upper 
left panel), heavy neutral scalar (upper right panel), pseudoscalar (lower left panel) and 
charged scalar (lower right panel) for the scalar masses and tan j3 of the benchmark scenario 
la defined in Table 4. 
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Table 11: Branching ratios as a function of cos(/3 — a) for the light neutral scalar (upper 
left panel), heavy neutral scalar (upper right panel), pseudoscalar (lower left panel) and 
charged scalar (lower right panel) for the scalar masses and tan j3 of the benchmark scenario 
lb defined in Table 4. 
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Table 12: Branching ratios as a function of cos(/3 — a) for the light neutral scalar (upper 
left panel), heavy neutral scalar (upper right panel), pseudoscalar (lower left panel) and 
charged scalar (lower right panel) for the scalar masses and tan j3 of the benchmark scenario 
2 defined in Table 5. 
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